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ABSTRACT 

An  experimental  investigation  was  conducted  to  determine 
the  dynamic  responses  of  stiffened  flat  plates  to  impulsive 
loading  by  underwater  shock,  waves.   Air-backed  flat  plates 
with  machined  external  stiffeners  and  clamped  boundary 
conditions  were  subjected  to  shock  loadings  from  TNT  charges 
detonated  underwater.   The  plates  were  instrumented  to 
measure  transient  strains  and  free  field  pressures.   Test 
results  were  compared  to  preshot  calculations  done  using  the 
USA-STAGS  code.   Particular  emphasis  was  placed  upon  the 
code 's ability  to  predict  stiffener  tripping  and  shear  at 
plate  boundaries. 
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I  .   INTRODUCTION 

A.   BACKGROUND  AND  HISTORY 

Sergei  Gorshkov,  Admiral  of  the  Soviet  Fleet,  has 
observed  that  the  end  of  World  War  Two  marked,  "... 
the  start  of  a  military-technical  revolution  which  in  scope 
and  depth  transcended  all  the  reforms  and  transmissions 
which  had  previously  occurred  in  the  armies  and  fleets  of 
the  world."  [Ref.  1] 

Much  has  been  made  of  the  tactical  and  geo-political 
implications  of  that  revolution.   Less  popularly  publicized, 
but  equally  radical  are  the  changes  that  the  developments  in 
technology  have  spawned  in  the  processes  of  ship  design  and 
construct  ion . 

Throughout  most  of  history,  war  at  sea  has  been  a 
process  of  trying  to  set  one's  opponent  afire,  board  him,  or 
hurl  chunks  of  stone  or  metal  at  him  across  relatively  short 
distances.   All  of  these  methods  had  the  distinguishing 
characteristic  of  being  essentially  disabling  tactics, 
unless  a  ship  was  burned  or  perforated  to  the  point  of 
sinking  (and  they  often  were  not),  damage  was  generally 
largely  superficial  and  the  ship  could  be  refitted  to  see 
service  again.   Witness  the  many  vessels  that  have  been 
fought  on  both  sides  of  a  conflict  and  the  tradition  of 
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expecting  naval  personnel  to  supplement  their  pay  with  prize 
money . 

Two  events  occurred  at  or  around  the  beginning  of  the 
twentieth  century  to  change  all  of  that:   the  invention  of 
high  explosives  and  the  advent  of  the  submarine. 

Initially,  the  submarine  was  primarily  a  commerce 
raider.   It  was  too  small  to  carry  a  prize  crew  and  too 
lightly  armed  to  effectively  combat  a  conventional  warship. 
During  World  War  One,  most  battles  between  submarines  and 
conventional  ships  were  fought  on  the  surface.   Submarines 
would  use  their  stealth  to  sneak-up  upon  unsuspecting  mer- 
chantmen, surface  and  then  attempt  to  sink  them  with  their 
deck  cannons.   Torpedoes  were  largely  unreliable  and  the 
submarines  themselves  were  at  risk  if  fast  armed  vessels 
were  present .   Convoying  the  merchants  with  small ,  agile 
men-of-war  was,  briefly,  an  effective  means  of  countering 
the  submarine  threat.   However,  as  torpedoes  became  more 
effective  so,  too,  did  the  submarine.   For  the  convoy  escort 
vessels  it  was  no  longer  sufficient  to  have  to  wait  for  a 
submarine  to  surface  in  order  to  attack  it.   An  effective 
standoff  weapon  was  needed. 

That  standoff  weapon  was  the  depth  charge.   First  intro- 
duced in  World  War  One,  the  depth  charge  came  of  age  in  the 
Second  World  War;  and  with  it  a  new  era  in  ship  design  was 
ushered  into  existence.   For  the  principle  of  this  new 
weapon  was  not  one  of  punching  a  hole  in  a  ship's  hull  and 
then  exploding  internally;  rather,  high  explosives  permitted 
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generation  of  intense  pressure  waves  in  the  sea  that  could 
break  a  vessel's  back  or  rupture  hull  fittings  whether  the 
vessel's  skin  was  punctured  or  not.   Simply  increasing  the 
thickness  of  plating  in  discrete  areas  of  vulnerability  was 
no  longer  sufficient.   Design  for  strength  and  resiliency 
under  unequal  pressure  loadings  was  required.   Moreover,' 
this  requirement  applied  to  attacker  as  well  as  attacked. 
Although  the  attacking  surface  vessel  had  the  advantage  of 
being  able  to  run  from  the  immediate  scene  of  the  explosion 
and  the  effects  were  not  compounded  by  as  great  a  pre- 
existing static  head  of  water,  these  weapons  were  not  direc- 
tional and  could  be  expected  to  damage  all  within  their 
sphere . 

The  issue  of  a  vessel's  ability  to  withstand  an  under- 
water explosion,  be  it  submarine  or  surface  ship,  became 
more  pressing  as  the  explosive  devices  became  increasingly 
powerful  and  sophisticated.   Not  only  did  depth  charges 
become  lethal  over  a  wider  area,  but  the  advancing  tech- 
nology permitted  mines,  too,  to  grow  from  essentially 
contact  weapons  to  devices  that  could  endanger  all  types  of 
shipping  from  considerable  depths.   The  issue  has  become 
particularly  crucial  with  the  development  of  nuclear  weapons 
which  are  so  powerful  that  no  device  need  contact  a  ship  or 
group  of  ships  to  cause  rampant  destruction. 

The  response  of  the  U.S.  Navy  was  to  require  that  ships 
and  their  critical  components  be  designed  and  built  to 
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withstand  specified  levels  of  shock  loading.   These  re- 
quirements included  (and  continue  to  include)  actual  shock 
testing  of  the  first  ship  of  a  class  and  many  crucial  systems 
on  the  component  level. 

During  the  immediate  post-World  War  Two  period,  many  of 
these  requirements  were  largely  tentative  and  experimental. 
By  the  early  1960 ' s,  however,  sufficient  experience  had  been 
gained  to  allow  design  requirements  to  be  set  forth 
[Ref.  2].   These  specifications  established  specific  guide- 
lines for  contractors  to  meet  based  upon  type  of  equipment, 
type  of  vessel  in  which  the  equipment  was  to  be  used,  and 
location  and  method  of  mounting  within  the  vessel.   As  with 
the  ships  themselves,  actual  shock  testing  of  the  item  was 
the  primary  means  of  establishing  compliance;  as  it  remains 
today . 

Performing  a  shock  test  upon  production  end  products-, 
particularly  when'  those  items  are  as  large  as  an  aircraft 
carrier  can  have  some  obvious  drawbacks:   it  is  costly,  it 
is  dangerous,  and  the  results  can  be  embarrassing. 
Moreover,  there  can  be  longterm  structural  damage  which 
is  neither  immediately  apparent  nor  correctable,  but  which 
does  not  by  itself  constitute  non-compliance  of  contract 
terms.   Moreover,  some  items  simply  do  not  lend  themselves 
to  shock  testing. 
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What  has  been  needed  is  a  method  of  analyzing  the  design 
itself  in  order  to  achieve  a  numerical  approximation  of  the 
end  product's  response  to  shock. 

This  need  has  coalesced  into  the  development  of  two 
computer  codes;  EPSA  (Elasto  Plastic  Shell  Analysis), 
produced  by  Weidlinger  Associates,  and  USA-STAGS  (Underwater 
Shock  Analysis  —  Structural  Analysis  of  General  Shells) 
produced  by  Lockheed  Missiles  and  Spacecraft,  Inc. 

B.   OBJECTIVES 

This  study  is  the  second  in  an  ongoing  series  sponsored 
by  the  Defense  Nuclear  Agency  (DNA)  into  the  effects  of 
underwater  explosions  upon  ship  hulls.   The  intent  of  these 
studies  has  been  to  test  the  applicability  of  the  two  previ- 
ously mentioned  structural  response,  underwater  shock  anal- 
ysis computer  codes  to  naval  design  and  analysis  requirements 
Under  the  umbrella  of  general  applicability,  focus  has 
been  placed  upon  two  defined  goals: 

1.  Use  of  the  codes  to  gain  insights  into  large  deflec- 
tion elastic-plastic  responses  of  a  submerged  struc- 
ture subjected  to  transient  acoustic  shock  loading; 
with  especial  emphasis  upon  stiffener  tripping 
phenomena . 

2.  Performance  of  underwater  explosion  testing  of  the 
structure  to  validate  the  results  predicted  by  the 
code  . 

To  achieve  these  goals,  flat  plate  models  geometrically 

similar  to  a  ship's  stiffened  hull  structure  have  been 

studied.   These  plates  were  air  backed  to  simulate  general 
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conditions  that  could  be  expected  in  a  ship's  plating;  and 
were  subjected  to  shock  loadings  that  produced  deformations 
well  into  the  plastic  regime  in  order  to  ensure  stiffener 
tripping.   Existing  information  concerning  underwater  explo- 
sion phenomena  was  employed  in  an  attempt  to  reduce  experi- 
mental uncertainties. 

The  first  study  in  this  series  was  that  conducted  by  LT 
Thomas  Rentz,  USN  in  fulfillment  of  thesis  requirements  at 
the  Naval  Postgraduate  School  [Ref.  3]  and  [Ref.  4].   In 
that  study  he  developed  the  following  experimental  strategy 
to  study  the  EPSA  code: 

An  attempt  to  validate  the  code  was  conducted  in  two 
phases:   Using  pressure  approximations  based  upon  empiri- 
cally determined  formulae,  the  code  was  first  used  as  a 
pre-test  predictor  of  test  results.   Once  the  underwater 
explosion  testing  had  been  conducted,  pressure  histories 
derived  from  that  testing  were  input  into  the  code.   The 
plate  strains  predicted  by  the  code  in  both  its  preshot  and 
post-shot  capacities  were  then  compared  to  actual  strain 
histories  from  testing.   Plate  deformations  and  boundary 
stresses  predicted  by  the  code  were  compared  to  observed 
physical  responses. 

In  this  study,  the  code  to  be  tested  is  USA-STAGS.   In 
an  effort  to  maintain  continuity  with  the  earlier  study  in 
the  series,  the  same  experimental  strategy  will  be  followed. 
The  test  geometry  and  flat  plate  model  will  also  be 
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retained.   Conclusions  regarding  the  ability  of  USA-STAGS  to 
predict  the  responses  of  a  simple  model  to  underwater  explo- 
sions will  be  drawn;  and,  perhaps  of  even  greater  impor- 
tance, guidelines  will  be  established  for  using  the 
USA-STAGS  code  and  conducting  underwater  explosion  testing 
for  future  studies. 
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II .   UNDERWATER  EXPLOSION  THEORY 

A.   THE  PRESSURE  PULSE 

Any  conversion  of  matter,  whether  chemical  or  nuclear, 
which  results  in  the  very  rapid  production  of  large  quanti- 
ties of  gas  at  very  high  temperature  can  be  categorized 
within  the  broad  species  "explosion."   Typically,  for  mili- 
tary high  explosives  such  as  TNT  the  pressure  and  tempera- 
ture within  the  gases  produced  by  an  explosion  are  on  the 
order  of  50,000  atmospheres  and  3000  degrees  Celsius 
[Ref.  5:   p.  3].   For  nuclear  explosions  the  initial  tempera- 
ture is  on  the  order  of  a  million  degrees  Celsius  and  the 
pressure,  which  is  dependent  upon  the  yield  is  similarly 
orders  of  magnitude  higher  than  that  for  a  conventional 
explosive  [Ref.  6]. 

Clearly,  temperatures  and  pressures  of  the  levels 
described  must  dissipate  high  levels  of  energy  through  any 
surrounding  medium.   In  the  case  of  water,  the  initial  mani- 
festation of  this  disturbance  is  an  intense  compression 
created  pressure  wave  propagating  radially  outward  from  the 
charge;  followed  nearly  instantaneously  by  a  displacement  of 
the  fluid.   In  the  immediate  vicinity  of  the  charge,  the 
velocity  of  propagation  of  the  pressure  wave  is  several 
times  the  speed  of  sound  (about  5000  feet  per  second  in 
water)  but  approaches  the  value  almost  immediately  and 
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therefore  may  be  apprehended  as  an  essentially  acoustic 
phenomenon.   As  the  water  begins  to  be  displaced,  the  gas 
bubble  created  by  the  explosion  expands;  thereby  decreasing 
its  internal  pressure  and  correspondingly  the  pressure  on 
the  constraining  fluid.   Viewed  on  a  pressure-time  curve  the 
initial  rise  in  pressure  is  so  nearly  vertical  as  to  appear 
discontinuous,  while  the  pressure  relief  is  characterized  by 
an  approximately  exponential  decay  (Fig.  2.1)  [Ref.  5:   pp. 
4-7]. 


Figure  2.1   Simple  Exponential  Approximation  of  Incident 
Shock  Wave 

The  actual  magnitudes  of  the  pressure  pulses  and  the  rate 

of  exponential  decay  have  been  found  through  experiment  to  be 
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functions  of  the  type  of  explosive,  the  charge  weight,  and 
the  standoff  distance.    The  equations  for  the  pressure  and 
decay  constant  were  originally  determined  by  Robert  H. 
Cole  and  were  presented  in  his  work,  Underwater  Explosions 
[Ref.  5].   These  have  since  been  simplified  to  use  a  series 
of  empirically  derived  coefficients  developed  by  R.  S. 
Price  [Ref .  7] . 

The  pressure  profile  can  be  expressed  by  the  equation 


t   -  t 

P(t )  =  P  exp  -^ 

o       o 


(eqn  2.1) 


where  P   is  the  initial  (greatest)  pressure  of  the  shock 
wave,  6  is  the  decay  constant  describing  the  exponential 
decay,  and  t  -  t1  is  the  time  elapsed  from  the  arrival  of 
the  shock  front . 

The  initial  pressure,  P  ,  and  the  decay  constant, 
9,  are  dependent  largely  upon  the  type  of  explosive  and 
the  weight  of  the  charge.   These  may  be  expressed 


P   =  K, 

o     1 


R 


(eqn  2.2) 


and , 


9  =  k2w 


1  /  ?\  A 
iff    \  2 

R 


(eqn  2.3) 


The  standoff  distance  is  defined  as  the  distance  from 
the  charge  to  the  nearest  point  on  the  target. 
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where  P   is  in  pounds  per  square  inch  and  9  in  milli- 
o 

seconds.   K  ,  K„  ,  A  ,  A„  are  the  explosive  dependent  values 
(Table  I).   W  and  R  are  the  charge  weight  (in  pounds)  and 
standoff  distance  (in  feet),  respectively. 

Table  I 
Explosive  Constants  I 

Explosive  TNT  HBX-1  PENTOLITE 

K  22505  22347.6  24589 

Ax  1.18  1.144  1.194 

K2  .058  .056  .052 

A„  -.185  -.247  -.257 


B.   SURFACE  CUTOFF  AND  BOTTOM  REFLECTION 

Shock  waves  propagate  radially  from  the  charge.   Thus, 
for  a  point  source  in  an  infinite  medium,  the  wave  could  be 
expected  to  travel  as  an  ever  expanding  sphere  until  damped 
by  the  fluid.   For  small  to  moderately  sized  charges  deto- 
nated well  below  the  surface  in  great  depths  of  water,  this 
is  a  satisfactory  analogy.   For  explosions  nearer  the 
surface  or  in  shoal  water,  however,  there  are  secondary 
effects  for  which  there  must  be  an  accounting. 

Vv'hen  a  shock  wave 'reaches  the  surface  of  the  water,  the 
adjoining  air  is  of  sufficiently  lesser  density  to  be  inca- 
pable of  supporting  the  shock  wave.  The  result  is  that  the 
effective  pressure  at  the  interface  must  be  zero  in 
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compliance  with  the  laws  of  continuity.   This  implies  that 
a  negat  ive  shock  wave  of  magnitude  equal  to  the  original 
incident  wave  must  be  propagated  downward  through  the 
medium  to  satisfy  equilibrium  conditions.   This  negative 
wave  will  travel  through  the  water  at  an  angle  to  the 
surface  that  is  equal  to  that  of  the  incident  wave;  much 
as  the  negative  reflection  image  from  a  mirror  continues  at 
an  angle  equal  to  the  angle  of  incidence. 

Remembering  that  these  waves  are  radiating  spherically 
outward  from  the  point  of  inception,  it  is  clear  that  both 
the  incident  wave  and  the  generated  surface  wave  will  even- 
tually reach  the  ship,  submarine,  or  other  target  object. 
It  is  likewise  clear  that,  with  its  longer  path,  the  surface 
wave  must  reach  the  target  at  some  time  later  than  the  inci- 
dent; which  travels  only  the  shortest,  or  standoff,  distance 
between  charge  and  target.   The  net  effect  is  additive. 
When  the  surface  wave  reaches  the  target,  the  initial  shock 
front  will  have  passed  and  the  incident  wave  will  be  at  some 
point  in  its  exponential  decay.   The  negative  surface  wave 
will  then  "chop  off"  the  tail  of  the  incident  wave.   If  the 
spatial  geometry  of  the  shot  is  such  that  the  magnitude  of 
the  surface  wave  is  greater  than  the  remaining  magnitude  of 
the  incident  wave,  a  region  of  negative  pressure  will  be 
created  (Fig.  2.2).   If  this  negative  is  less  than 
the  pressure  required  to  keep  the  fluid  in  a  liquid  state, 
the  fluid  will  "flash"  into  vapor.   This  chopping  effect  is 
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referred  to  as  surface  cutoff,  and  the  negative  pressure  as 
bulk  cavitation. 
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Figure  2.2   Exponential  Approximation  with  Surface  Cutoff 

The  direct  effect  of  bulk  cavitation  is,  by  itself,  not 
particularly  harmful.   The  structure  simply  experiences 
essentially  no  pressure  while  it  is  in  the  cavitated  region 
When,  however,  the  static  head  of  water  above  the  region 
overcomes  the  cavitation  forces,  the  region  closes  suddenly 
to  generate  a  destructive  reloading  upon  the  structure. 

To  further  complicate  matters,  if  the  shot  occurs  in 
shoal  waters  or  at  depths  near  the  bottom  a  second  boundary 
is  also  present.   Unlike  the  surface,  in  this  case  the 
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incident  wave  is  impinging  upon  a  medium  of  greater  density 
than  the  fluid.   The  reflected  wave  will  therefore  be 
positive.   The  actual  magnitude  of  this  wave  depends  upon 
the  physical  characteristics  of  the  bottom  material.   The 
effects  of  this  bottom  reflected  wave  are  also  additive 
and  will  generally  reach  the  target  at  some  time  after 
the  initial  shock  front  has  passed;  although  the  actual 
time  of  arrival  of  this  secondary  front  will  once  again 
be  dependent  upon  the  bottom  conditions  and  the  speed  at 
which  the  wave  travels  through  the  bottom  medium.   Addi- 
tionally, in  shoal  waters,  that  same  bottom  reflected  wave 
will  reach  the  surface;  creating  its  own  surface  cutoff. 

The  net  effect  upon  the  target  is  thus  one  of  an  inci- 
dent shock  wave;  which,  if  the  charge  and  target  are  suffi- 
ciently near  the  surface,  will  be  cutoff  by  the  negative 
surface  wave.   This,  in  turn,  is  followed  by  the  bottom 
reflected  wave  which  also  experiences  an  exponential  decay, 
and  which  may  also  be  cutoff  by  its  own  generated  surface 
wave  ( Fig .  2.3). 

Analytically,  the  easiest  way  to  model  the  effects  of 
incident,  surface,  and  bottom  waves  upon  the  target  is  as 
has  been  represented  in  Figure  2.4,  with  three  separate 
images . 

The  incident  and  surface  images  are  relatively  easy  to 
depict.   The  incident  image,  M. ,  of  course  represents  the 
charge  itself  and  occupies  the  same  location  in  the  spatial 
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Figure  2.3   Complex  Shock  Wave  Pattern 


25 


M. 


*0-_ 


R 


SURFACE 


R> 


BOTTOM 


M 


b  O 


Figure  2.4   Three  Image  Presentation  of  Shock  Waves 
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geometry.   The  surface  image,  M  ,  is  of  the  same  magnitude 
as  the  incident  image  but  of  an  opposite  sign.   It  lies  on  a 
line  drawn  vertically  through  the  charge  at  a  distance  above 
the  surface  equal  to  the  depth  of  the  charge  below. 

The  surface  cutoff  pressure  may  therefore  be  computed 

pc  =  -Ki(r— )  (eqn  2-4) 

where  R   is  now  the  distance  from  the  surface  image  to 
s  ° 

the  target.   The  time  of  cutoff  after  the  passage  of  the 

initial  shock  front  is  thus  the  time  difference  for  the  wave 

to  cover  the  greater  distance,  or 

R   -  R. 
s     i 

t   =  (eqn  2.5) 

c       c  H 

where  c  is  the  speed  of  sound  in  the  fluid.   The  cutoff 
pressure  will  then  decay  in  the  same  exponential  manner  as 
the  incident  pressure. 

As  may  be  conceived  from  the  preceding  discussion  of  the 
bottom  reflection;  the  location  of  the  bottom  image  is  not 
so  straightforward.   It  may  be  viewed  as  initiating  along 
the  same  vertical  line  as  do  the  incident  and  surface 
images,  and  at  a  depth  below  the  fluid-bottom  interface 
equal  to  the  distance  of  the  charge  above.   Here,  however, 
the  similarities  to  the  surface  image  end.   The  bottom  image 
cannot  be  considered  a  point  charge  because  of  the  tendency 
of  the  pressure  wave  to  "smear"  across  the  bottom.   The 
magnitude  and  speed  of  the  wave,  too,  are  dependent  upon 
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bottom  characteristics.   A  muddy  bottom,  for  instance,  will 
tend  to  reflect  a  diffuse,  low  magnitude  wave  travelling  at 
or  about  the  speed  of  sound  in  water.   A  rock  bottom  on  the 
other  hand  will  produce  sharp,  high  magnitude  pulses  which 
will  tend  to  radiate  over  a  distance  due  to  the  higher  speed 
of  travel  of  the  shock  wave  in  the  bottom  medium  than  in  the 
fluid.'  Due  to  this  dependence  upon  bottom  conditions, 
determination  of  bottom  reflection  effects  do  not  lend  them- 
selves to  general  predictive  studies  and  are  often  disre- 
garded, although  they  can  be  significant. 

C.   BUBBLE  EFFECT 

The  foregoing  discussion  of  the  development  of  the  pres- 
sure wave  and  its  impact  upon  the  target  assumes  that  a 
single,  clean  wavefront  is  produced  by  the  explosion.   In 
point  of  fact,  such  is  rarely  the  case.   When  the  charge 
explodes,  it  first  creates  a  high  intensity  shock,  front 
which  is  the  initial  incident  shock  wave.   It  also  produces 
a  bubble  of  very  hot  gasses  at  very  high  pressure. 

The  bubble  expands  until  it  reaches  pressure  equilibrium 
with  the  surrounding  water.   When  it  does  so,  a  portion  of 
the  energy  within  the  bubble  is  released  as  a  new  shock 
front.   This  release  of  energy  then  permits  the  momentum  of 
the  water  displaced  to  assert  itself  and  the  bubble  is 
collapsed,  compressing  the  gases  within.   The  displaced 
water  can  be  conceived  as  attempting  to  regain  its  original 
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pre-bubble  condition.   As  the  water  rushes  in,  it  compresses 
the  gases  beyond  the  equilibrium  condition  and  the  bubble 
again  expands,  again  releasing  a  portion  of  its  energy  as  a 
shock  wave,  albeit  a  smaller  one  than  the  preceding  (Fig. 
2.5).   This  process  has  been  known  to  produce  as  many  as  ten 
significant  pulses  [Ref.  5:   pp.  8,  9],  which,  though  weaker 
than  the  initial  shock  front,  can  still  have  a  pronounced 
effect  upon  a  structure  which  has  already  been  loaded  by  the 
earlier  pulses . 


FIRST 
MAXIMUM 


SECONO 
MAXIMUM 


FIRST    PERIOD 


— t- — SECOND    PERIOD 


Figure    2.5      Multiple    Expansions    of    the    Hot    Gas    Bubble 
[Ref.    5 :       p .    7 ] 
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The  time  and  radius  of  the  bubble  at  the  first  pulse 
has,  again,  been  determined  empirically,  the  maximum  radius 
as 

f        w1/3'   \ 

A     =  Kj T-T^\  (eqn  2.6) 

3V(D  +  33)1/3j 

at  time, 

(      w1/3     A 

t  =  K   =-^  (eqn  2.7) 

H(D  +  33)°/6y 

where  A     is  expressed  in  feet  and  t  in  seconds.   W  is 
max 

once  more  the  weight  of  the  charge,  K„  and  K.  are  explosive 
dependent  constants  (Table  II),  and  D  is  the  depth  of  the 
charge  in  feet . 


Table  II 
Explosive  Constants  II 


Explosive  HBX-1  TNT  PENTOLITE 

K  14.14  12.67  12.88 

KA  4.761  4.268  4.339 


A  second  and  distinct  phenomenon  is  bubble  migration;  so 
called  because  of  the  tendency  of  the  bubble  to  migrate 
toward  any  nearby  mass.   When  the  bubble  contacts  the  mass, 
it  collapses  and  begins  a  series  of  rapid  pulsations.   The 
net  effect  is  similar  to  that  of  the  already  described 
bubble  pulse;  a  series  of  reloadings  upon  the  structure. 
The  rapidity  of  the  pulsations  of  the  collapsed  bubble  and 
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their  increased  local izat ion  tend,  however,  to  magnify  their 
severity  . 

D.   APPLICATIONS  OF  THEORY  TO  EXPERIMENT  DESIGN 

Developing  a  basic  understanding  of  the  primary 
phenomena  involved  in  underwater  explosion  theory  is 
crucial  to  developing  a  successful  underwater  test  program 
for  many  reasons.   Two  are  of  particular  interest  here, 
however : 

First,  and  most  obviously,  equations  2.1  through  2.3 
allow  initial  estimates  of  the  peak  pressure  and  pressure 
history  to  which  a  test  structure  is  to  be  exposed.   These 
are  then  input  to  the  computer  code  so  that  preshot  approxi- 
mations of  the  structural  responses  may  be  made;  thereby 
creating  guidelines  for  calibrating  gages,  determining 
charge  weights,  etc. 

Second,  the  design  of  the  test  itself  is  a  direct  conse- 
quence of  the  phenomena  being  observed.   As  is  readily 
apparent  from  the  earlier  sections  of  this  chapter,  the 
loadings  imposed  upon  a  structure  by  an  underwater  explosion 
can  be  extremely  complex.   The  incident  wave  may  be  a 
simple  exponential  decay  dependent  only  upon  charge  and 
standoff,  but  the  many  secondary  effects  are  subject  to 
variables  which  cannot  be  readily  controlled  experimentally. 

The  contributions  of  the  secondary  effects  are  signifi- 
cant and  must  eventually  be  studied.  However,  during  these 
early  phases  of  investigation  into  underwater  explosion 
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effects  a  determined  effort  has  been  made  to  simplify  the 
problem  as  much  as  possible. 

It  has  therefore  been  required  that : 

1.  The  model  to  be  studied  be  as  basic  as  is  consis- 
tent with  stated  objectives. 

2.  Study  be  limited  only  to  the  incident  shock  wave. 
These  two  requirements  are  interrelated  and  are  tied 

intrinsically  to  the  nature  of  underwater  explosion 
phenomena.   As  has  been  seen  earlier,  the  incident  shock  may 
be  approximated  as  a  plane  wave  emanating  from  the  charge 
and  striking  the  panel.   In  the  absence  of  secondary 
effects,  the  responses  being  studied  here  may  thus  be 
limited  to  those  of  the  test  panel  alone.   If,  on  the  other 
hand,  secondary  effects  are  included,  more  than  a  single 
point  source  is  involved  and  the  interactions  of  any  support 
or  backing  structures  with  the  fluid  medium  must  also  be 
considered . 

Fortunately,  for  the  purposes  of  this  study,  surface 
cutoff  and  bottom  reflection  are  later  time  responses  which 
can  be  neglected  if  study  is  confined  to  the  first  few 
milliseconds  after  detonation.   Bubble  and  cavitation 
effects  can  be  reduced  greatly  or  eliminated  through  careful 
attention  to  the  test  geometry. 

If  an  explosion  occurs  sufficiently  near  the  surface 
that  the  initial  bubble  radius  is  greater  than  the  charge 
depth,  the  hot  gases  forming  the  bubble  will  vent  into  the 
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atmosphere  and  be  dissipated.  This  venting  is  instant aneous 
and,  if  occurring  at  the  proper  depths,  prevents  any  further 
bubble  pulsation  or  migration. 

Bulk  cavitation,  too,  is  dependent  upon  the  test  geom- 
etry.  The  region  of  cavitation  tends  to  expand  horizontally 
away  from  the  line  drawn  between  source  and  surface  images 
(Fig.  2.6).   The  actual  shape  and  dimensions  of  the  cavita- 
tion region  will  vary  with  charge  size,  type,  and  depth,  but 
the  area  at  depths  below  the  charge  can  generally  be 
expected  to  remain  free  of  cavitation  effects. 
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Region  of  Bulk  Cavitation 

(From  Weidlinger  Associates  CUE  Code.) 
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Bulk  cavitation  and  bubble  effects  can  therefore  be 
reduced  significantly  or  eliminated  entirely  if  the  test 
geometry  is  such  that  the  charge  is  located  directly  above 
the  plate  at  a  depth  to  cause  bubble  venting. 

In  designing  the  test  platform  [Ref.  3:   pp.  19-28]  a 
series  of  iterations  were  run  using  the  EPSA  code  to  deter- 
mine a  charge  weight  and  standoff  distance  that  would 

produce  maximum  test  panel  deflections  of  approximately  four 

2 

plate  thicknesses.    From  these  runs  it  was  determined  that 

the  ideal  charge  weight  and  standoff  were  8  pounds  of  TNT  at 
9  feet. 

The  maximum  bubble  radius  is  expressed  as  a  function  of 
charge  weight  and  depth  (Eqn.  2.6).   Since  it  has  been  found 

that  bubble  venting  occurs  when  the  ratio  of  charge  depth  to 

3 

maximum  radius  is  less  than  0.75,   choosing  a  charge  depth 

is  a  fairly  straightforward  procedure  once  a  determination 
of  the  desired  ratio  has  been  made. 


2 

This  value  was  recommended  by  Weidlinger  Associates  as 

one  which  experience  indicated  would  ensure  stiffener 
tripping . 

3 

Any  ratio  below  0.75  will  cause  venting.   Venting  by 

its  very  nature  will  duct  a  portion  of  t he  .explosion  energy 
into  the  atmosphere,  thereby  decreasing  the  maximum  pressure 
of  the  shock  wave.   As  a  rule  of  thumb,  this  decrease  in 
usually  considered  negligible  until  the  ratio  drops  below 
0.50. 
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By  selecting  a  depth  to  radius  ratio  of  0.50  it  was 
readily  determined  that  a  charge  depth  of  four  feet  was 
appropriate.   Applying  this  to  Eqn .  2.6  yields 

(            R1/3 
A     =  (12.67) — Q 

\(4  +  33)1/2 


=  7.60  ft 
4  ft 


7.60  ft 


0.526 


which  assures  venting. 

From  the  foregoing  it  can  be  seen  that  a  test  configura- 
tion that  will  have  the  desired  consequences  of  reducing 
bubble  and  bulk  cavitation  effects  can  be  designed.   For  an 
eight  pound  TNT  charge,  this  platform  will  have  the  charge 
centered  vertically  at  a  standoff  distance  of  nine  feet 
above  the  plate.   The  charge  itself  will  be  located  at  a 
depth  of  four  feet  below  that  surface.   The  actual  design  of 
the  test  platform  will  be  discussed  more  fully  in  Chapter  V. 
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III.   THE  FLAT  PLATE  MODEL 

As  stated  in  the  introduction,  the  objective  of  this 
study  has  been  to  test  the  applicability  of  a  computer  code, 
USA-STAGS,  to  predicting  the  effects  of  an  underwater  explo- 
sion upon  ships'  hulls.   To  do  this,  it  is  clearly  necessary 

that  some  basis  for  modelling  simple  sections  of  a  typical 

4 
ship's  grillage   be  established. 

The  scaling  laws  used  here  and  in  the  earlier  work  of 
this  series  [Ref.  3]  were  those  developed  by  Dr.  Raymond  P. 
Daddazio  of  Weidlinger  Associates,  Inc. 

In  his  work  [Ref.  8],  Dr.  Daddazio  has  established  two 
dimensionless  parameters,  3  and  \,    which  he  refers  to  as 
the  "plate  slenderness  ratio,"  and  the  "longitudinal  stif- 
fener  slenderness  ratio,"  respectively.   These  are  defined 


o  b   _y  (eqn  3.1) 

p     t   E 


and 


x  -  &  TT  (ec>n  3-2) 


4 
The  ship's  hull  or  superstructural  plating  and  associ- 
ated stiffeners. 
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where 

a   : 
y 

E 
a 
b 
t 
k 


=  material  yield  strength 

=  Young's  modulus 

=  spacing  between  transverse  frames 

=  spacing  between  longitudinal  stiffeners 

=  plate  thickness 

=  radius  of  gyration  of  the  longitudinal  stiffener 

acting  with  an  assumed  effective  width  of  plating 

b   where , 
e       ' 


K     k  2     1 

6      6  "  ^2 


Through  sampling  a  variety  of  ship  types,  Dr.  Daddazio 
found  that  typical  values  for  ships'  grillage  fall  into  the 
range : 

1-.0  <  6  <  4.5 

and  , 

0.15  <  X    <    0.9 

The  intent  of  this  series  of  studies  has  been  to  use  one 
basic  model  but  with  a  variety  of  stiffener  types  so  that 
the  responses  of  different  stiffener  geometries  could  be 
studied. 

The  model  established  in  Ref.  3  was  a  3/16  inch  thick 
test  panel  18  inches  long  by  12  inches  wide  milled  into  the 
center  of  a  two  inch  thick  6061-T6  aluminum  blank  measuring 


37 


33  inches  by  27  inches.   Two  siiffeners  with  rectangular 
cross  sections,  3/16  inch  thick  by  one  inch  deep  were  mounted 
transversely  across  the  test  panel  in  such  a  manner  as  to  be 
symmetric  to  both  the  transverse  and  longitudinal  axes  (Fig. 
3.1).   The  entire  assembly  was  machined  as  a  single  unit  to 
avoid  the  inconsistencies  of  welding.   A  radius  was  cut  at 
all  corners  and  the  whole  polished  to  avoid  stress  concen- 
tration points. 

To  apply  Equation  3.1  and  Equation  3.2  the  two  milled 
stiffeners  were  considered  to  be  longitudinal  stiffening 
members  and  the  longer  boundaries  of  the  panel  to  be 
transverse  frame  members.   Thus,  the  dimensions  of  the  plate 
are,  a  =  12  inches,  b  =  6  inches,  and  t  =  3/16  inches.   Using 
the  standard  reference  values  for  6061-T6  aluminum,  it  was 
found  that,  for  this  original  plate:   [Ref.  3] 

6  =  2.02 

and 

\    =    0.9 

which  fall  within  the  acceptable  ranges. 

The  initial  plate  used  in  this  study  was  identical  to 
that  used  in  Ref.  3  with  the  exception  that  the  stiffener 
cross  section  used  was  a  tee  rather  than  a  rectangle  (Fig. 
3.2)  . 

Since  the  spacing  of  the  stiffeners  in  this  tee- 
stiffened  plate  was  identical  to  that  in  the  original  plate, 
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Figure  3.1 


Test  Panel  with  Transverse  Rectangular  Stiffeners 
[Ref.  3:   p.  13] 
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Figure  3.2   Test  Panel  with  Transverse  Tee  Stiffeners 
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clearly  the  3  values  remained  invariant .   Given  ihe 
changed  geometry  of  the  stiffener  ihe  effective  radius  of 
gyration  of  the  plate /st if fener  assembly  (k)  could  be 
expected  to  change,  and  with  it  the  value  of  X. 

Before  that  value  X  could  be  determined,  however,  the 
dimensions  of  the  stiffener  itself  first  had  to  be 
established . 

In  an  effort  to  ensure  that  stiffener  responses  were  a 

function  of  stiffener  design  alone,  the  plastic  moduli,  Z  , 

of  the  two  stiffener  cross  sections  were  set  as  close' as 

possible  to  equal.   For  the  original,  rectangular-stiffened 

plate,  Z   equaled  0.046875  in  . 
p 

By  a  process  of  iteration,  a  tee-st if fener  with  flange 
dimensions  23/32  inch  by  1/8  inch,  and  web  dimensions  25/32 
inch  by  1/8  inch  was  fixed  upon  as  giving  an  elastic  modulus, 
Z   =  0.046170  in3,  which  is  within  1.5%  of  that  for  the 
rectangular  stiffener  (Fig.  3.3). 

Now,  applying  these  dimensions  to  the  grillage  param- 
eter, X,  we  find  that  the  neutral  axis  of  the  effect  ive 
stiffener  (Fig.  3.4)  is: 


y  = 


EA.y. 

r  i 


0.2220  in 


The  moment  of  inertia  of  the  effective  stiffener  is 


L„  =  ZI  .  +  ZA  .  d  .  ' 
ZZ      i      ii 

=  0.092562  in 


4 
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Plastic  Section  Modulus: 

Z   =  Sum  of  the  moments  about  the 
neutral  axis  of  the  cross 
sections  above  and  below  the 
axis 


Z     =   (b.h.)(h,  +  hh.) 

p  112  1 


h,    = 


1/8    in.    =   0.125000    in 
23 


h„    = 


32 

25 

32 


in.    =   0.718750    in. 
-   0.607747 


-   0.173503    in. 


h      =   0.607747    in. 

b0    =   b.    =  ^  in.    =   0.  125000    in 
2  3         o 


Z      =    (0.718750) (0. 125) (0.173503   +  H   ~   0.125) 
P 

+  H    (0.125) (0.173503)2 


+  k    (0.  125)  (0.607747) 
0.046170    in3 
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Figure    3.3      Plastic    Modulus    of    Tee-Stiffened    Plate 
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Figure  3.4   Effective  Cross  Section  of  Tee-St if f end  Plate 
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which  gives  a  radius  of  curvature: 


k  = 


A 


-  [0.092562  in4/1.0256  in"]2 


=  0.300419  in 


Thus,  for  the  tee-stiffened  plate, 


X    = 


12 


40000 


30041971   10  x  106 


=  0.804144 

which  is  well  within  the  parameter  range  required  for 
modelling  ships'  grillage. 

The  second  requirement  levied  upon  these  test  plates  was 
the  ability  to  demonstrate  stiffener  tripping. 

Stiffener  tripping  has  been  identified  as  the  lateral- 
torsional  instability  of  the  stiffener  as  it  becomes 
suddenly  unstable  and  fails  under  load.   It  is  identified  by 
a  characteristic  warping  and  buckling  of  the  stiffener.   It 
is  significant  that  the  continued  resistance  of  a  stiffened 
structure  to  deformation  under  loading  is  dramatically 
reduced  after  the  stiffeners  have  tripped.   In  earlier 
studies  of  impulsive  loadings  upon  ships'  grillage,  tripping 
has  been  identified  as  "...  a  primary  ductile  failure  mode 
for  ship  structure."  [Ref.  9:   p.  2]. 

Previous  observations  have  indicated  that  tripping  is 
most  readily  initiated  when  an  axially  compressive  loading 


44 


is  applied  to  the  top  fibers  of  the  stiffener.   Recognizing 
that  a  flat  plate  will  tend  to  deform  away  from  the  source 
of  an  explosion,  the  plates  in  this  series  of  studies  have 

been  oriented  with  the  stiffener  sides  outward  into  the 

5 
fluid  medium. 

As  stated  earlier  in  this  chapter,  the  original  intent 

in  this  series  of  studies  was  to  maintain  an  essentially 

fixed  plate  model,  varying  only  the  stiffener  geometries. 

Unfortunately,  the  results  of  the  first  two  explosion  tests 

(on  the  original  rectangular-stiffened  plate  [Ref.  3],  and 

on  the  tee-stiffened  plate  here)  were  not  satisfactory  in 

view  of  the  stated  objectives.   The  two  transversely  mounted 

stiffeners  were  too  short  and  the  entire  assembly  too  rigid 

to  ef feet ively. demonstrate  the  desired  stiffener  tripping. 

The  shock  wave  essentially  "punched"  through  the  panel 

material,  shearing  it  at  the  fixed  boundaries  and  at  the 

6 
base  of  the  stiffeners   (Fig.  3.5). 

The  results  of  the  two  initial  tests  indicated  that  a 

design  change  was  in  order.   The  third  plate  of  the  series 

(the  second  of  this  study)  was  a  simplified  model  to  permit 


5 
This  is  hardly  the  expected  orientation  on  a  surface 

vessel,  but  is  common  on  many  submarine  pressure  hulls. 

6 
Although  not  easy  to  see  in  Figure  3 . d ,  not  all  of  the 

damage  to  the  tee-stiffened  plate  was  done  by  the  incident 

shock  wave.   Apparently,  when  the  bottom  reflection  hit  the 

plate  and  its  backing  structure,  it  was  forced  upward  into 

its  own  supporting  cables.   Much  of  the  damage  to  the  stif- 

finers  in  particular  seems  to  have  been  caused  by  the 

actions  of  that  cable. 
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Figure  3.5   Explosion  Damage  to  Tee-Stiffened  Plate 
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closer  investigation  of  the  shock  effects  themselves, 
without  attempting  either  to  retain  a  set  plastic  modulus 
or  to  approximate  ships'  grillage. 

The  result  of  this  simplification  was  a  thickening  of 
the  base  panel  to  J  inch  to  increase  resistance  to  shear. 
The  stiffener  was  restored  to  the  original  rectangular 

configuration  but  the  thickness  was  increased  to  i  inch 

7 
and  the  depth  to  If  inches.    Most  important,  in  this 

attempt,  a  single,  longitudinally  oriented  stiffener  was 
used  rather  than -the  two  transverse  stiff eners  in  an  effort 
to  provide  greater  flexibility  along  the  length  of  the 
stiffener  (Fig.  3.6). 

Although  the  charge  used  to  test  the  rectangular  stif- 
fened plate  was  small  (|  pound  of  TNT),  and  the  correspond- 
ing plate  deflections  were  far  too  slight  for  any  firm 
conclusions  to  be  drawn,  it  can  be  seen  in  Figures  3.7  and 
3.8  that  tripping  action  appears  to  have  begun  to  manifest 
itself.   The  single,  longitudinally  oriented  stiffener  seems 
to  be  the  superior  model;  and,  while  further  testing  is 
required  at  larger  charge  weights ,  its  employment  is  tenta- 
tively recommended  for  future  studies. 


This  decision  was  made  somewhat  arbitrarily.   Dimen- 
sions were  chosen,  in  part,  to  reduce  the  cost  and  simplify 
machining  the  plate.   It  is  interesting  to  note,  however, 
that,  with  3  =  1.52  and  \  =  1.1,  this  plate  is  not  far 
outside  the  acceptable  grillage  range. 
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Figure  3.6   Test  Panel  with  Longitudinal  Rectangular  Stiffener 
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IV.   THE  USA-STAGS  CODE 

A.   GENERAL 

Structural  Analysis  of  General  Shells  (STAGS)  is  the 
finite  element  code  first  developed  by  Lockheed  Missiles  and 
Space  Company,  Inc.  during  the  1960 ' s  in  response  to  an 
initial  tasking  by  the  Navy  and  the  Air  Force  plus  an  inter- 
nally perceived  need  for  a  general  structural  analysis 
program. 

The  current  version  of  STAGS  is  STAGS  C-l ,  a  dynamic 
analysis  code  with  broad  applications  to  the  aerospace 
industry  which  has  also  found  application  in  ship  design  and 
engineering.   A  new  version,  STAGS  2000,  is  now  under  develop- 
ment which  will  have  an  entirely  changed  data  manager 
concept  from  that  used  in  the  C-l  version.   This  change 
should  make  the  new  code  both  easier  to  use  and  far  more 
economical  of  the  user's  time.   STAGS  2000  will  appear 
initially  as  a  static  analysis  code,  but  will  eventually 
exceed  STAGS  C-l's  capabilities. 

STAGS  C-l  is  supported  by  Computer  Software  Management 
and  Information  Center  (COSMIC)  under  the  sponsorship  of  the 
National  Aeronautics  and  Space  Administration  (NASA). 

The  USA-STAGS  code  is  an  extension  of  the  basic  STAGS 
routine  to  permit  evaluation  of  impulsive  loadings  upon 
design  structures  in  a  fluid  medium.   As  such,  Underwater 
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Shock  Analysis  (USA)  is  a  mull imodule  routine  which  is 
inserted  as  an  intermediary  and  adjunct  to  the  program  pre- 
analyzer  (STAGS1)  and  analyzer  (STAGS2).   The  STAGS1 ,  STAGS2 
combination  can  stand  alone  for  some  types  of  general  struc- 
tural analysis.   As  a  result  of  its  add-on  nature,  the  USA 
code  cannot  be  used  in  a  stand  alone  mode,  but  it  is  not 
STAGS  dependent.   Versions  have  been  created  for  use  with 

Q 

other  finite  element  structural  analysis  codes. 

B .   THEORY 

The  characteristics  which  make  USA-STAGS  unique  lie  in 
its  applicability  to  a  wide  variety  of  problems  and,  partic- 
ularly, in  the  individual  techniques  used  to  solve  those 
problems.   The  broad,  underlying  theory  behind  those  solu- 
tions is  in  itself  common  to  most  finite  element  codes. 

Briefly,  the  technique  used  by  codes  of  this  genre  is  to 
solve  the  differential  equations  of  motion  for  a  structure 
in  order  to  determine  its  responses  to  its  environment 
[Ref.  10:  pp.  2-1  through  2-14].   By  then  applying  the 
constitutive  equations  (Hooke's  law  generalized  to  include 
anisotropic  materials),  the  stresses  in  the  governing  equa- 
tions of  motion  may  be  expressed  in  terms  of  strains.   These 
strains  may  in  turn  be  expressed  as  the  first  partials  of 
the  displacements  with  respect  to  the  spatial  coordinates 


As  of  this  writing,  USA  has  been  linked  to  NASTRAN , 
GENSAM,  SPAR,  and  DIAL. 
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through  application  of  the  kinematic  relationships.   An 
equivalency  can  be  drawn  between  the  motions  of  the  struc- 
ture between  any  two  instants  in  time  t~  and  t1  and  the 

extrema  of  the  energy  functional  /  1(U  -  K)dt  with  respect 

0 
to  all  admissible  variations  of  the  strain  energy  U  and  the 

kinetic  energy  K  which  satisfy  the  geometric  boundary  condi- 
tions and  continuity  conditions.   The  equilibrium  equations 
can  then  be  derived  from  the  first  variation  of  the 
energy  functional  through  partial  integration  with  respect 
to  time  and  the  spatial  variables. 

The  variational  approach  provides  a  particularly  effec- 
tive means  for  analyzing  continuous  systems  because  of  the 
way  in  which  it  generates  and  accounts  for  natural  boundary 
conditions  [Ref.  11].   The  disadvantage  of  the  method  is 
that  it  is  limited  to  conservative  systems.   However,  by 
applying  the  principle  of  virtual  work  STAGS  is  able  to 
greatly  extend  the  classes  of  problems  to  which  it  can  be 
applied . 

In  the  second  case;  where  inelastic  deformations  such  as 
plastic  strains  or  thermal  expansions  are  encountered,  STAGS 
treats  them  as  "pseudo-loads."   In  this  manner,  by  applying 
these  pseudo-load  terms  to  the  force  side  of  the  equation, 
an  initial  estimate  of  the  inelastic  deformation  values  can 
be  made  and,  if  convergence  is  obtained,  the  nonconservat ive 
problem  can  be  solved  by  iteration  upon  a  series  of 
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conservative  problems  [Ref.  10:   pp.  2-7,  2-8].   The  range 
of  cases  to  which  STAGS  may  now  be  applied  is  shown  in 
Figure  4.1. 

After  being  reduced  to  functions  of  displacement  only, 
the  governing  equations  of  motion  for  a  discretized  system 
may  be  expressed  in  the  form: 

M  X  +D  (y)  +B  (•<)  +K  (x)  =  f 
s~     s  ~     s-     s~ 

where  x  is  the  vector  of  displacement  components,  M   is 

the  structural  mass  matrix,  f  is  the  vector  of  external 

forces,  and  K   is  the  (generally  nonlinear)  structural 

stiffness  operator.   The  operators  B   and  D   include 

s       s 

forces  that  are  functions  of  structural  deformation  and 
deformation  velocity,  respectively. 

Equation  4.1  is  the  matrix  equation  that  STAGS  must 
solve  to  find  structural  displacements.   For  fluid/struc- 
tural interactions  the  generalized  force  vector  becomes  a 
much  more  complex  entity  [Ref.  12:   pp.  2-1  through  2-12]. 

For  submerged  structures  excited  by  an  acoustic  wave, 
the  exciting  forces  are  given  by: 


f  =  -GAf(PJ  +  P  )  (eqn  4.2) 

where  PT  and  PQ  are  the  incident  and  scattered  nodal  pres- 
sure  vectors  acting  upon  the  surface  of  the  structure  at  the 
fluid-structural  interface.   Af  is  the  diagonal  area  matrix 
associated  with  elements  in  the  fluid  mesh,  and  G  is  the 
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matrix  which  relates  the  structural  and  fluid  surface 
forces.   It  can  be  seen  from  the  development  of  the  pressure 
characteristics  of  an  underwater  shock  wave  in  Chapter  1, 
that  for  an  underwater  explosion  the  incident  pressure 
vector  PT  is  known.   The  scattered  pressure  vector  PQ  on  the 
other  hand  is  dependent  upon  the  configuration  of  the  struc- 
tural surface  and  is  therefore  not  only  not  immediately 
known  but  is  also  changing  as  the  structure  deforms. 

It  is  to  solve  this  computational  dilemma  and  thereby 
solve  the  combined  equation  of  motion: 


u   x  +  Ds(x)  +  Bs<x)  +  Ks(x)  =  _GAf(?s  +  ?s} 


(eqn  4.3) 


that  the  USA  code  was  created. 

The  USA  code  finds  PQ  by  utilizing  the  Doubly  Asymptotic 
Approximation  (DAA)  and  a  staggered  solution  procedure. 

The  DAA  is  an  approximate  relation  which  approaches 
exactness  at  both  the  high-frequency  (early-time)  and  low- 
frequency  (late-time)  limits  to  the  passage  of  the  shock 
front.   It  may  be  written: 


MfPs  +  pcAfPs  =  pcMfv_ 


(eqn  4.4) 


where  vQ  is  the  vector  of  scattered-wave  fluid-particle 
velocities  normal  to  the  structural  surface  at  the  fluid- 
structural  interface.   Mf  is  the  fluid  mass  matrix  for  the 
fluid  mesh  at  the  interface,  and  p  and  c  are  the  density  and 
sound  velocity  of  the  fluid,  respectively. 
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When  fluid  and  structure  are  excited  by  an  acoustic 
wave,  the  incident  and  scattered  fluid  particle  velocities 
(v   and  vQ)  may  be  related  to  the  structural  response  by  the 
kinematic  compat ability  relation: 

T* 
G  x  =  VT  +  vq  (eqn  4.5) 

Solving  for  v„  and  substituting  the  resultant  into  equation 
4 . 4  yields 

M,F>   +  pcA„PQ  =  pcM„(G  x  -  v,)  (eqn  4.6) 

Equations  4.3  and  4.6  are  now  each  expressed  in  terms  of  the 

single  unknown  PQ  and  a  solution  may  be  obtained  by  a  stag- 

gered  iteration  process  between  the  two  equations  for  each 
t ime  step . 

C.   ORGANIZATION 

As  has  already  been  stated,  USA-STAGS  is  a  modular  code. 
Functionally,  it  may  be  divided  into  five  basic  components. 

1.  The  structural  preprocessor 

2.  The  fluid  mass  preprocessor 

3.  The  augmented  matrix  preprocessor 

4.  The  time  integration  processor 

5.  The  data  postprocessor 

Each  of  these  components  is  run  individually  in  sequence 
utilizing  the  outputs  from  the  previous  modules. 
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1 .   The  Structural  Preprocessor 

The  structural  preprocessor  is  the  unit  that 
assembles  the  structural  mass  and  stiffness  matrices, 
creates  the  finite  element  mesh,  and  otherwise  provides  an 
encoded  description  of  the  structure's  internal  and  ex- 
ternal degrees  of  freedom  [Ref.  12:   pp.  3-3,  3-4]. 

For  USA-STAGS,  the  structural  preprocessor  is  the 
STAGS1  portion  of  the  STAGS  code  itself. 

Input  for  STAGS  is  made  to  each  of  ten  functional 
subdivisions  which  correspond  to  the  types  of  input  'that  go 
into  each.   These  subdivisions  are: 

-  Summary  and  Control  Parameters 

-  Computational  Strategy  Parameters 

-  Data  Tables 

-  Geometry  and  Material 

-  Discretization 

-  Discrete  Stiffeners* 

-  Boundary  Conditions 

-  Loads* 

-  Output  Control 

-  Finite  Element  Units* 

(Subdirectories  designated  with  *  are  not  required  with  all 
models . ) 

Input  to  each  subdirectory  is  in  the  form  of  a 
series  of  data  records  consisting  of  real  and  integer 
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numerical  data  fields.   This  will  be  discussed  in  greater 
detail  in  the  following  section. 

2 .  The  Fluid  Mass  Preprocessor 

The  fluid  mass  preprocessor  (FLUMAS)  constructs  the 
fluid  mass  matrix  for  a  structure  submerged  in  a  fluid. 
Inputs  include  [Ref.  12:   pp.  3-4,  3-5]  the  global  coordi- 
nates for  the  structural  geometry  (input  from  the  structural 
preprocessor),  plus 

-  Fluid  wet-surface  mesh 

-  Element  definitions 

-  Material  properties  of  the  fluid 

-  Location  of  free  surface 

-  Model  symmetries 

3 .  The  Augmented  Matrix  Preprocessor 

The  augmented  matrix  preprocessor  (AUGMAT)  receives 
the  output  from  the  structural  and  fluid  mass  processors  to 
create  the  matrices  required  for  solution  of  the  augmented 
equations  4.3  and  4.6. 

4 .  The  Time  Integration  Processor 

The  functional  heart  of  the  USA-STAGS  system,  the 
time  integrator  processor  (TIMINT)  links  back  to  and  supple- 
ments the  structural  processor  STAGS2  to  implement  the  stag- 
gered solution  technique;  creating  the  solution  data  files 
over  the  time  range  specified.   Input  includes  the  data 
passed  from  AUGMAT  plus 
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-  Spatial  geometry  of  charge  and  target 

-  Initial  incident  pressure  and  pressure  wave  decay 
characteristics 

-  Time  integration  information 

-  Restart  data 

TIMINT  includes  an  internal  postprocessing  capability. 
Generally,  however,  this  function  is  left  to  the  following 
module . 

5 .   The  Data  Postprocessor 

USA-STAGS  Provides  a  selection  of  three  possible 
postprocessors:   POSTP ,  POSTPR,  and  STAPL .   Their  uses  will 
be  discussed  in  the  following  section. 

D.   USING  THE  CODE 

As  has  already  been  stated,  the  STAGS  code  is  extensive 
in  its  concept;  attempting  to  cover  in  a  single  code  as 
many  approaches  to  problem  solution  as  can  be  linked  in 
a  single  contextual  framework. 

It  is  beyond  the  scope  of  this  work  to  cover  all 
possible  aspects  of  the  code;  indeed  within  the  stated  capa- 
bilities of  STAGS  this  approach  may  seem  peripheral,  being 
limited  as  it  is  to  those  aspects  clearly  related  to  under- 
water shock  analysis  of  flat  plate  models. 

STAGS  is  designed  to  be  used  for  any  one  of  six  types  of 
analysis  [Ref.  13:   p.  2-1]. 

-  Linear  stress  analysis 

-  Geometrically  nonlinear  elastic  stress  analysis 
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-  Inelastic  stress  analysis,  geometrically  linear  or 
nonl inear 

-  Bifurcation  buckling  analysis  with  linear  or  geometri- 
cally nonlinear  stress  state  (elastic) 

-  Small  vibration  analysis  based  on  linear  or  geometri- 
cally nonlinear  stress  state  (elastic) 

-  Transient  response  analysis,  linear  or  geometrically 
nonlinear,  elastic  or  inelastic 

In  this  study,  only  two  of  those  capabilities  were  used: 
(a)  Small  vibrations  (in  a  stress  free  state),  and  (b) 
Geometrically  nonlinear  transient  response  analysis.   Of 
these,  the  second  alone  is  truly  applicable  to  USA-STAGS. 
The  small  vibration  analysis  was  used  as  a  check  on  the  code 
and  input  data.   Since  the  actual  inputs  for  the  two  anal- 
yses vary  but  slightly,  both  shall  be  addressed  together. 

Before  that  can  be  done,  however,  the  codes  must  somehow 
be  brought  into  some  semblance  of  conformity;  for  they  are 
not  normally  provided  as  a  single  unit  . 

Both  codes  are  provided  in  tape  form:   STAGS  from      , 
COSMIC,  and  USA  directly  from  Lockheed.   These  are  designed 
to  be  loaded  onto  the  VAX  11/780  computer  using  the  standard 
MOUNT  MF:   command.9 

It  is  suggested  that  ,  at  least  during  the  early  stages 
of  processing,  each  tape  be  assigned  to  a  separate  direc- 
tory, for  many  of  the  source  files  will  be  culled  as  the 


9 
See  Appendix  A. 
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executable  files  ( * . EXE  )    are  created.   The  source  files 
supplied  with  each  code  are  specific  to  the  code  until  the 
final  linkage,  and  their  sheer  number  can  cause  confusion  if 
not  approached  methodically. 

The  STAGS  tape  as  provided  contains  all  the  necessary 
ingredients  for  the  functioning  codes  except  the  executables 
themselves.   These  are  left  to  be  linked  on  the  individual 
computer.   The  link  commands  for  creation  of  these  executa- 
bles are  provided  with  the  tape  on  the  INFO. VAX  file,  which 
also  contains  a  full  listing  of  the  tape's  contents.   The 
commands  that  are  of  immediate  interest  are  LSTAGS1.COM, 
LSTAGS2.COM,  LP0STP.COM,  and  LSTAPL.COM.   These  should  be 
copied  into  the  file  and  executed  in  the  normal  manner 
for  a  program  run  from  the  terminal.     The  result  will  be 
the  executable  files  for  the  respective  module. 

Executables  are  provided  with  the  USA  tape;  however, 

since  relinking  may  be  required  to  adapt  the  code  to  local 

12 
filenames,  link  command  files  are  also  provided.     The 


On  the  VAX,  the  asterisk  represents  a  "wildcard"  which 
may  be  substituted  for  any  filename,  filetype,  or  version 
number  (the  complete  listing  for  a  file  is:   filename, 
filetype,  version  number)  whenever  a  general  category  is 
desired  or  the  listing  is  to  be  repeated. 

This  is  done  by  simply  typing  "@  command  filename"  and 
hitting  the  return  key. 

12 

The  USA  tape  also  includes  another  particularly  nice 

feature,  the  complete  user's  manual  for  each  of  the  modules. 
This  ensures  that  the  most  current  version  is  available; 
and,  more  important,  the  version  of  the  manual  that  is 
pertinent  to  that  tape.   A  vital  accompaniment  since  the 
tapes  are  continually  being  updated. 
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command  files  to  be  used  in  constructing  the  USA  executables 
are  LINKFLU.COM,  LINKAUG.COM,  LINKUSAS.COM,  and  LINKP0S.COM. 
It  is  the  LINKUSAS.COM  file  which  actually  performs  the 
marriage  between  the  codes;  allowing  the  integration/ 
structural  analysis  functions  of  STAGS  to  be  coupled  with 
the  USA  fluid  processor. 

Once  all  of  the  executable  files  have  been  created  and 
tested  in  a  full  run  of  the  combined  code;  all  files  but 

the  executables  may  be  deleted,  thus  realizing  a  considerable 

,        13 
increase  in  computational  space. 

The  basic  strategy  that  will  be  followed  here  shall  be 
to  address  USA-STAGS  as  a  single  code  broken  into  five 
different  modules  (see  preceding  section);  each  of  which  can 
be  run  independently  in  stepping  stone  fashion  once  neces- 
sary input  from  earlier  modules  has  been  obtained. 

The  first  module  to  be  run  is  STAGS1 .   The  actual  input 
for  this  includes  all  information  needed  for  both  STAGS1  and 
STAGS2  execution  in  a  single  entry  file. 

The  first  five  data  records:   A-l ,  B-l ,  B-2 ,  B-3,  and 
C-l  are  common  to  all  STAGS  runs.   The  inputs  as  set  forth 
in  the  COSMIC  STAGS  users  guide  [Ref.  13:   pp.  3-1  through 
3-11]  are  largely  self-explanatory  but  a  few  points  should 
be  noted. 


13 

Total  space  required  for  the  executables  alone  is  less 

than  4000  blocks;  compared  to  the  almost  13000  blocks  for 

the  executables  plus  source  codes. 
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First,  IPOSTO,  IPOST1  ,  and  IPOST2  are  applicable  only  to 
the  STAGS  postprocessors  POSTP ,  and  STAPL .   If  the  user  is 
interested  only  in  the  nodal  displacements  resulting  from 
underwater  shock  these  may  be  derived  directly  from  either 
of  the  USA  modules  TIMINT  or  POSTPR .   Thus  IPOSTO,  IP0ST1  , 
and  IPOST2  may  be  set  to  zero  thereby  suppressing  creation 
of  model  and  solution  data  files  (*.MOD  and  FOR022.DAT)  and 
realizing  some  savings  in  computer  space.   It  is  suggested 
that  this  be  done  in  any  case  during  early  model  development 
and  while  testing  code  inputs  (as  with  the  small  vibrations 
analysis).   If,  however,  strains  are  required;  these  are 
computed  only  within  the  STAGS  postprocessors  and  IPOSTO 
must  be  set  to  1.   Considerable  amounts  of  computer  time  and 
space  can  be  used  in  processing  a  large  model;  time  which 
would  be  wasted  if  this  is  not  done.   Similarly,  IP0ST1  and 
IP0ST2  should  be  set  to  a  non-zero  value  if  displacement  and 
stress  plots  are  desired  at  the  IPOSTlst  and  IP0ST2nd  load 
steps . 

Second,  it  should  also  be  noted  well  that  IFLU  must  be 
set  to  1  if  underwater  shock  analysis  using  USA  is  to  be 
undertaken.   This  will  cause  the  FOR003.DAT  and  FOR004.DAT 
files  required  for  USA  to  be  created.   These  files  are  not 


14 

This  is  required  only  when  using  STAPL.   POSTP  is 

capable  of  calculating  displacements  and  stresses 

internally . 
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unduly  lengthy,  but  for  the  sake  of  economy  should  be 
suppressed  for  non-USA  runs . 

The  C-l  data  record  is  meaningless  for  small  vibration 

analysis  and  transient  response  analysis  but  must  be 

15 
included.   A  value  of  0.  should  be  entered. 

The  only  true  variation  in  data  record  sequence  between 
small  vibration  analysis  and  transient  response  analysis 
comes  in  the  next  group  of  cards.   The  proper  sequencing  for 
small  vibrations  is  D-2  ,  D-3,  F-l;  while  for  transient  anal- 
ysis it  is  D-l,  E-l ,  E-2 ,  F-l.   These  are  the  Strategy 
Parameters  and  the  difference  lies  in  the  need  for  integra- 
tion over  time  in  the  transient  model.   The  D-l  record  sets 
tolerances,  limits  the  amount  of  CPU  time,  allows  restarts, 

etc.-;  while  the  E  records  control  the  actual  time  integra- 

16 
tion  itself. 
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This  must  be  a  real  number.   STAGS  is  often  blind  to 

data  not  submitted  in  the  proper  format.   It  will  not  record 

an  error,  but  will  simply  pass  over  the  improperly  input 

data  until  finding  some  in  the  format  it  seeks.   This  can 

lead  to  sometimes  dangerous  distortions  of  the  input 

records . 

16 

At  the  time  that  this  study  was  initiated  it  was 

understood  that  implicit  integration  using  K.  C.  Park's 
formula  ( IMPL  =  0,  METHOD  =  4  on  the  E-2  record)  was  the 
only  form  of  integration  that  should  be  used  with  USA.   As 
of  this  writing,  however,  Lockheed  informs  that  USA  may  now 
also  be  used  with  explicit  integration.   It  has  been 
suggested  that  this  might  provide  a  faster  and  more  econom- 
ical solution.   Due  to  time  limitations  that  has  not  been 
tested  here. 
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The  D-2  and  D-3  records  are  specific  to  modal  analysis 
and  deal  with  the  number  and  range  of  eigenvalues  to  be 
calculated . 

Data  records  F  through  P  pertain  to  the  model  itself. 
The  records  F  through  L  provide  information  general  to  the 
entire  model  and  are  listed  in  normal  alphabetical  sequence 
once  through  only,  although  specific  cards  may  be  repeated 
within  the  sequence.   M  through  P,  plus  Q  and  R  are  specific 
to  the  individual  shell  and  beam  components  and  are  repeated 
in  sequence  once  for  each  component  (see  examples  in  the 
appendix )  . 

The  F-l  and  F-2  records  provide  grid  point  and  stiffener 
summaries,  respectively.   The  information  in  F-l  controls 
generation  of  the  finite  element  mesh.   Inputs  are  simply 
the  number  of  nodal  rows  and  columns  for  each  shell.   Actual 
choice  of  which  shell  direction  shall  consist  of  rows  and 
which  columns  is  unimportant  as  long  as  that  choice  remains 
consistent  for  the  local  X  and  Y  coordinates  across  the 
entire  model  (rows  are  lines  of  consistent  X,  columns  are 
lines  of  consistent  Y).   Entry  is  made  sequentially;  number 
of  rows ( I ) ,  number  of  columns(I)  for  each  shell  in  order  of 
input.   The  mesh  created  forms  the  master  elements  which 
will  then  be  fitted  into  the  structural  geometry  by  the 
information  contained  in  the  M  records. 

The  F-2  record  is  a  summary  of  shell  stiffeners.   Use  of 
this  record  could  be  extremely  useful  in  determining 
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generalized  responses  on  large  scale  models,  but  was  incon- 
sistent with  the  purposes  of  this  study  because  of  the 
inherent  assumption  that  stiffeners  do  not  deform. 

Record  G-l  establishes  the  compatibility  conditions  for 
the  various  shells.   The  instructions  listed  in  the 
COSMIC-STAGS  user's  manual  are  straightforward,  but  refer- 
ence to  Appendix  B  may  prove  helpful  for  establishing  proper 
linkages.   It  will  be  noted  that  in  each  case,  for  a  flat 
plate,  side  1  is  that  nearest  and  parallel  to  the  global  Y 
axis  while  side  4  has  a  similar  correspondence  to  the  global 
X  axis . 

Records  I,  J,  K,  and  L  constitute  the  Data  Table  section 
of  STAGS.   Again  the  user's  manual  is  self-explanatory  and 
will  not  be  explicated  except  to  note  that  the  Material 
Record,  1-1;  the  Material  Property  Record,  1-2;  and  the 
Plasticity  Record  1-3  (if  used);  are  listed  in  sequence 
once  for  each  of  the  materials  enumerated  in  the  B-3  record. 

If  plastic  behavior  is  expected  in  the  model  and  the 
White-Besseling  theory  of  plasticity  is  chosen  ( IPLST 
options  0  or  1  on  the  1-1  record)  especial  note  should  be 
made  of  the  user's  manual's  warnings  concerning  the  slope  of 
the  line  connecting  points  submitted  for  the  discretized 
stress-strain  curve  in  the  1-3  record.   The  change  in  slope 
at  each  point  must  be  negative.   Another  point  not  mentioned 
in  the  user's  manual,  but  which  should  definitely  be  noted, 
is  that  the  first  stress-strain  data  entry  on  the  1-3  record 
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must  be  consistent  with  the  elastic  modulus  entered  on  the 
1-2  record.   That  is,  the  slope  of  a  line  drawn  from  the 
origin  to  that  first  data  entry  should  be  the  elastic 
modulus.   Also,  while  specifying  plastic  behavior 
may  actually  decrease  the  run  time  for  the  TIMINT  module  by 
reducing  the  number  of  iterations  required  for  each  time 
step,  it  can  increase  the  required  computer  space  by  more 
than  100%. 17 

If  plasticity  is  to  be  included  in  the  model,  there  are 
also  some  changes  that  should  be  made  in  the  K  records. 

STAGS  determines  the  plastic  responses  of  the  shell  by 
integrating  through  the  thickness.   The  input  NLIP  in  the 
K-2  record  establishes  the  number  of  integration  points 
through  each  material  layer.   Since  Simpson's  rule  is  used, 
this  number  must  be  odd .   NLIP  may  lie  between  0  and  9 
(inclusive),  but  some  care  must  be  exercised  in  the  choice, 
since  the  larger  numbers  will  add  to  the  length  of  the 


17 

A  normal  run  for  a  model  with  about  300  nodes  seems  to 

be  on  the  order  of  about  two  to  two-and-a-half  hours  of  CPU 

time  for  strictly  elastic  behavior.   If  the  elastic  range  is 

exceeded  but  plastic  behavior  is  not  specified,  these  runs 

can  easily  continue  for  15  to  20  hours  as  the  code  attempts 

repeated  iterations  for  each  time  step  (the  user  can  control 

this  somewhat  by  specifying  the  maximum  allowable  time  for 

any  iteration),  and  the  full  number  of  specified  time  steps 

may  not  be  attainable.   Inclusion  of  plasticity  will  again 

reduce  the  run  time  to  roughly  three  to  four  hours.   Output 

of  stresses  beyond  the  yieM  point  of  the  material,  however, 

requires  much  larger  solution  data  (POR022.DAT)  files  for 

both  the  TIMINT  and  POSTP  runs.   Typically,  these  have  been 
on  the  order  of  45,000  and  35,000  blocks  respectively  for 
the  pressure  levels  used  here. 
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TIMINT  run;  which,  as  has  already  been  mentioned,  is  not 
insignificant  for  the  plastic  case.   The  default  value  for 
plasticity  is  3.   This  may  be  acceptable  for  models  where 
the  membrane  stress  governs,  but  will  not  be  so  if  there  is 
marked  bending  stress. 

Another  change  that  should  be  made  for  runs  including 
plastic  responses  is  to  enter  a  value  of  1  or  2  for  the  LSO 
entry  of  the  K-2  record.   If  this  is  not  done  the  USA-STAGS 
code  will  successfully  complete  an  entire  run,  but  the  post- 

1  Pi 

processor,  POSTP ,  will  not  print  stress  or  strain  data. 

The  sequence  of  records  M  through  R  defines  the  proper- 
ties and  loads  on  the  shell  units.   As  is  noted  most  explic- 
itly in  the  user's  manual,  these  records  are  read 
sequentially  for  each  unit.   All  data  in  the  sequence  will 
be  read  before  the  M-l  record  for  the  next  unit  is  read. 

The  individual  shell  geometries  and  the  defining  parame- 
ters by  which  they  are  to  be  connected  to  other  shells  (that 
is  the  means  by  which  they  are  to  be  oriented  to  the  global 
geometry)  are  specified  in  record  M-l.   IGLOBE  option  0  is 
useful  if  the  local  coordinates  of  the  unit  already  corre- 
spond precisely  to  the  global  coordinates.   Option  1  is 


1 8 

Actually,  this  may  not  be  entirely  true.   In  order  to 

get  any  useful  response  whatsoever  out  of  POSTP,  it  has  been 

necessary  to  "trick"  the  code  into  producing  the  desired 

stresses  and  strains  by  running  the  module  with  an 

FOR002.DAT  file  created  by  a  STAGS1  run  for  the  elastic 

case.   This  has  been  true  whether  or  not  LSO  has  been 

flagged,  with  no  noticable  difference  in  the  results. 
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similarly  intended  to  be  a  labor  saving  device.   A  caveat 
should  be  observed,  however,  that  option  1  should  only  be 
used  if  the  side  connecting  two  shells  is  uniquely  defined 
(i.e.,  no  straight  lines).   It  is  strongly  recommended  thai 
options  2  or  3 ,  which  allow  more  user  control  by  specifying 
boundaries  through  specific  global  reference  points, be  used 
in  preference  to  option  1. 

An  option  not  listed  in  the  user's  manual,  but  which 
could  conceivably  be  useful  for  construction  of  cylindrical 
models  is  option  4.   This  allows  the  user  to  specify  the 
location  of  a  unit  origin  in  the  global  coordinate  system 
and  then  define  the  boundary  in  terms  of  one  translation  and 
one  rotation  about  the  reference  point. 

It  should  also  be  noted  that  STAGS  will  permit  introduc- 
tion of  shell  geometries  not  included  in  the  eleven  which 
are  indigenous  to  the  program.   This  is  done  through  a  user 
written  subroutine  which  may  be  specified  by  choosing  ISHELL 
option  1  on  the  M-l  record. 

If  the  ISHELL  =  1  option  is  chosen,  the  shell  geometric 
properties  normally  input  in  record  M-2  will  be  input  in  the 
user-written  subroutine  LAME  as  stated  above.   Otherwise, 
shell  geometries  are  chosen  from  those  provided  by  the  code. 
Property  values  for  these  geometries  should  be  input  as 
specified  in  the  figures  of  Appendix  B.   If  plastic  deforma- 
tions are  expected,  particular  attention  must  be  paid  to 
inputs  to  the  M-5  record.   First,  since  Gaussian  quadrature 
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is  used  to  integrate  across  the  mesh,  STAGS  will  give  erro- 
neous results  (or  none  at  all)  if  that  mesh  is  not  located 
on  the  midplane  of  the  shell  thickness.   ECZ  must  corre- 
spondingly be  set  to  0,  and  IPLAST  set  equal  to  1.   No 
further  inputs  need  be  made  on  this  record  unless  one  of  the 
shell  units  is  located  on  a  plane  of  symmetry.   If  that  is 
the  case,  STAGS  will  not  compute  lateral  deformations 
without  further  information  about  how  the  shell  is  to 
respond.   This  information  is  provided  in  the  form  of  minute 
randomly  generated  displacements  or  "random  imperfections" 
of  the  shell  nodes  from  the  midplane.   To  trigger  input  of 
these  imperfections,  I RAMP  should  be  set  equal  to  1,  and  the 
appropriate  entries  made  on  the  M-6  and  M-7A  records. 

The  user's  manual's  instructions  and  flow  patterns  for 
records  N  through  R  are  clear  and  it  is  believed  that,  with 
possible  references  to  the  input  files  in  Appendix  C,  the  user 
should  have  no  difficulty  completing  the  STAGS  input  deck. 

As  stated  in  the  section  on  organization  of  the  code, 
USA  is  divided  into  four  modules.   Each  of  these  may  be  run 
independently  once  the  needed  data  files  have  been  produced 
by  the  previous  module.   Except  for  the  required  access  to 
those  generated  data  files,  each  module  is  run  entirely  in 
stand-alone  fashion.   Since  modular  input  information  is 
brief  and  error  statements  are  generally  good,  this  permits 
debugging  without  affecting  earlier  inputs  and  generates  a 
high  degree  of  confidence  in  the  modular  output. 
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It  has  already  been  mentioned  that  setting  the  IFLU 
switch  equal  to  1  on  the  B-l  record  of  the  STAGS1  run  will 
generate  two  data  files  FOR003.DAT  and  FOR004.DAT  which  are 
specific  to  USA.   These  are  produced  in  addition  to  the 
FOR002.DAT  and  output  files  which  are  normally  generated  by 
a  STAGS1  run. 

Each  of  the  three  data  files  FOR002.DAT,  FOR003.DAT, 
and  FOR004.DAT  must  be  renamed  after  the  STAGS1  run  for  use 
in  USA.   The  first,  FOR002.DAT,  should  simply  be  renamed  for 
the  file  subdirectory  containing  the  USA  files  (e.g., 
(*.USA)FOR002.DAT) .   FOR003.DAT  and  FOR004.DAT  must  be 
renamed  to  correspond  to  the  GRDNAM  and  MASNAM  files  speci- 
fied in  the  FLUMAS  and  AUGMAT  modules,  respectively.   This 
may  be  done  by  entering  a  simple  RENAME  command  on  the  VAX. 
The  user  will  then  be  queried  regarding  the  name  from  which 

and  to  which  they  are  being  changed. 

The  GRDNAM  file  (the  renamed  FOR003.DAT  file)  contains 
the  global  coordinates  of  the  structural  grid  points  gener- 
ated for  the  model  by  STAGS1 .   The  FLUMAS  module  then  uses 
this  information  to  construct  a  fluid  mass  matrix  (FLUNAM) 
as  well  as  generate  the  fluid  mesh  data  (GEONAM)  and  trans- 
formation coordinates  that  relate  the  structural  and  fluid 
degrees  of  freedom  on  the  wetted  surface.   The  FLUNAM  and 
GEONAM  data  files  will  then  be  used  in  the  AUGMAT  module  to 
create  a  final  augmented  mass  matrix. 
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The  module  AUGMAT  is  the  component  of  USA  which  assem- 
bles the  data  from  the  FLU.MAS  module  and  from  STAGS1  into 
the  specific  constants  and  arrays  required  for  the  staggered 
solution  process.   The  input  deck  for  this  module  is  very 
short  and  experience  indicates  that  problems  with  the  module 
will  generally  be  caused  by  input  errors  to  either  STAGS1  or 
FLUMAS . 

USA's  time  integrator  module,  TIMINT,  links  the  results 
of  the  AUGMAT  module  which  are  stored  in  the  data  file 
PRENAM  back  to  STAGS2  and  then  unites  with  STAGS2  to  conduct 
a  step-by-step  numerical  time  integration  of  the  governing 
equations  for  submerged  structures  exposed  to  shock  waves. 
The  time  integration  information  input  to  TIMINT  will  over- 
write earlier  entries  to  the  STAGS  deck.   It  is  suggested, 
however,  that  these  entries  should  be  made  consistent. 

Input  to  TIMINT  includes  not  only  time  integration 
information  but  also  the  spatial  coordinates  for  the  charge 
and  for  the  point  on  the  structure  closest  to  the  charge 
(the  standoff).   These  distances  are  made  in  reference  to 
the  origin  of  the  global  structural  coordinate  as  defined  in 
the  STAGS  input.   One  point  of  caution:   The  surface-cutoff 
provision  in  FLUMAS  assumes  a  horizontal  separation  between 
charge  and  target.   As  the  charge  approaches  a  point  verti- 
cally above  the  target,  some  of  the  computed  values  will  go 
to  infinity,  causing  the  module  to  fail.   The  surface-cutoff 
function  does  not  automatically  turn  itself  off  and  the 
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FLUMAS  input  rausi  be  changed  accordingly  if  a  vertical 
aspect  is  required. 

Provision  has  also  been  made  in  TIMINT  to  print  nodal 
displacements  without  recourse  to  a  post-processor.   It  is 
recommended  that  this  option  not  be  exercised.   A  run  for 
any  of  the  other  USA-STAGS  modules  can  be  made  in  a  matter 
of  minutes.   Running  TIMINT  can  easily  consume  several  hours 
of  CPU  time  for  a  moderately  sized  model.   It  is  therefore 
not  time  or  cost  effective  to  run  TIMINT  each  time  a  new  set 
of  displacements  is  needed. 

The  solution  data  file,  POSNAM ,  output  by  TIMINT  is 
specific  to  the  USA  post-processor,  POSTPR.   If  nodal 
displacements  only  are  desired  this  is  by  far  the  easiest 
and  most  economical  route  to  follow.   It  is  not,  however, 
the  only  option.   USA-STAGS  provides  a  selection  of  three 
possible  post-processors.   As  mentioned,  POSTPR  is  internal 
to  USA  and  may  be  used  to  provide  either  data  listings  or 
graphical  displays  of  nodal  displacements.   The  other  two, 

POSTP  and  STAPL ,  are  a  part  of  the  STAGS  package.   A  second 

19 
solution  data  file  in  STAGS  format,   FOR022.DAT,  is  created 


19It  is  actually  created  by  STAGS2  during  the  TIMINT 


run 
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and  may  be  renamed  to  the  STAGS  subdirectory  for  use  in 

20 

either  of  these. 

Both  POSTP  and  STAPL  have  been  designed  to  provide  not 
only  displacements  but  also  stresses  and  strains;  POSTP  in 
tabular  form  and  STAPL  in  graphical  on  a  CALCOMP  plotter. 
At  the  time  of  this  writing,  the  Naval  Postgraduate  School 
did  not  have  a  graphical  capability  on  VAX /VMS  so  neither 

the  POSTPR-graphical  nor  STAPL  options  were  run  at  the 

21 
Naval  Postgraduate  School. 

In  accordance  with  common  usage  on  the  VAX,  the  sequence 
of  commands  required  to  execute  any  module  have  been  entered 
into  command  files  ( * . COM ) .   Examples  are  given  in  Appendix 
A. 

Once  the  appropriate  file  names  have  been  entered  into 
the  command  file  for  the  module  to  be  run,  the  command  file 
may  be  executed  by  entering  "@  command  filename"  for  runs 
from  the  terminal;  or  "SUBMIT  filename"  for  batch  runs.   A 
file  definition  statement  (SET  DEF  [*.*])  is  required  for 
batch  runs  but  is  not  necessary  for  runs  from  the  terminal. 


20It  should  be  noted  that  the  TIMINT  entry,  NSAVER,  is 
the  frequency  with  which  responses  are  saved  for  the  USA 
postprocessor  only.  The  frequency  with  which  data  is  saved 
for  the  POSTP  and  STAPL  postprocessors  is  actually  tied  to 
the  restart  data  in  NRESET .  If  stress  and  strain  informa- 
tion is  required  at  each  timestep,  NRESET  should  therefore 
be  set  to  1. 

21 

POSTPR  and  STAPL  plots  shown  in  later  sections  of  this 

work  were  done  by  Lockheed. 
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E.   DISCRETIZING  THE  MODEL 

As  its  name  implies,  the  primary  building  unit  in  the 

STAGS  code  is  the  shell  element.   These  elements  may  be 

constructed  into  shapes  which  are  warped  through  almost  any 

conceivable  configuration  in  three-space  and  beam  elements 

may  be  added,  but  the  basic  unit  remains  two-dimensional. 

The  flat  plate  models  constructed  for  this  study,  then, 
were  ideally  suited  for  discretization  using  such  elements. 

Each  plate  with  its  associated  stiffeners  was  capable  of 
being  divided  into  several  flat  shell  units  composed  of 
square  or  nearly  square  elements.   Additionally,  the  even 
spacing  of  the  stiffeners  allowed  advantage  to  be  taken  of 
the  plate  symmetries  to  reduce  the  overall  size  of  the  model 
and  corresponding  length  of  the  computer  run. 

The  transverse  tee-stiffened  plate  established  symmetry 
in  two  directions  permitting  a  quarter  model  to  be  used. 
That  quarter  model  was  then  divided  into  five  shells;  two 
for  the  base  plate,  and  three  for  the  stiffener.   These 
shells  were  numbered  beginning  with  the  two  shells 
comprising  the  base  upward  through  the  stiffener  (Fig. 
4.2).   This  was  done  to  permit  changes  to  be  made  to  the 
code  input  by  simply  adding  or  deleting  shells.   As  long  as 
the  basic  plate  dimensions  and  symmetries  remained  the  same, 
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extensive  rewriting  could  be  avoided  throughout  the 

22 

USA-STAGS  modules. 

The  STAGS  code  will  generate  the  finite  element  mesh. 
The  user  need  only  specify  the  number  of  nodal  rows  and 
columns.   In  the  case  of  the  tee  stiffened  plate,  an  element 
0.6  inches  by  0.6  inches  was  chosen  in  order  to  keep  the 
mesh  as  fine  as  possible  and  yet  be  able  to  place  the  active 

filaments  of  the  strain  gages  used  during  the  underwater 

23 

explosion  testing  entirely  within  the  element.     This  corre- 
sponded to  11  nodal  rows  and  6  nodal  columns  in  shell 
unit  1,  11  rows  and  11  columns  in  shell  unit  2,  and  11 
rows  and  3  columns  in  each  of  the  three  stiffener  shells 
(Figure  4.3). 

Since  the  stiffener  used  on  the  second  test  plate  was 
rectangular,  three  shell  units  sufficed  to  model  it  (Fig. 
4.4).   At  the  time  that  this  plate  stiffener  was  discre- 
tized,  the  limit  set  upon  the  number  fluid  mesh  elements  by 


the  FLUMAS  module  was  171. 


24 


The  increased  surface  area 


22 

Unfortunately,  it  was  not  possible  to  carry  through 

with  this  plan  since  the  plate  with  the  longitudinal  rectan- 
gular stiffener  completely  changed  material  dimensions  and 
permitted  only  one  axis  of  symmetry. 

23 

STAGS  computes  stresses  and  strains  at  the  geometric 

centroids  of  the  elements.   Strain  gages  should  therefore  be 
centered  on  the  element . 

24 

It  turns  out  that  this  is  relatively  easy  to  change, 

but  it  involves  technical  considerations  that  are  beyond  the 

scope  of  this  work.   Questions  about  such  modification 

should  be  directed  to  Lockheed.   The  limit  on  most  accounts 

has  recently  been  increased  to  400. 
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necessitated  by  the  use  of  a  halt-plate  model  instead  of  a 
quarter-plate  set  the  total  number  of  fluid  elements  beyond 
that  number  when  the  0.6  inch  element  was  used.   Rectangular 
elements  1.2000  inches  by  1.2857  inches  were  chosen  for 
the  base  shells  in  order  to  allow  strain  gage  placement  at 
the  center  of  each  shell.   Elements  for  the  stiffener  shell 

ere  1.2857  inches  by  0.6250  inch.   This  corresponded  to  6 
nodal  rows  and  8  nodal  columns  for  shell  units  1  and  2 , 
d  3  rows  and  8  columns  for  shell  unit  3  (Fig.  4.5). 


w 


an 


Figure  7T~  Half  Model  of  Rectangular-Stiffened  Plate 
For  the  final  model  of  both  the  tee-stiffened  and 
rectangular-stiffened  plates,  the  finite  element  mesh  was 
placed  at  the  mid-planes  of  the  shell  units.   This  was  not 
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originally  done.   In  an  attempt  to  make  the  USA-STAGS 
output  correspond  directly  to  the  strain  gage  outputs  from 
underwater  testing,  the  plates  were  modelled  with  the  meshes 
placed  on  the  surfaces  of  the  plates  which  corresponded  to 
the  surface  of  strain  gage  mounting.   The  half-thickness  of 
the  plate  was  then  input  as  the  ECZ  value  on  the  M-5  record, 
thereby  defining  the  true  mid-plane  for  the  code.   This 
procedure  worked  extremely  well  for  elastic  runs.   As  stated 
in  the  previous  section,  however,  the  Gaussian  quadrature 
that  STAGS  uses  to  calculate  the  plastic  responses  across 

the  mesh  surface  is  incapable  of  handling  the  integrations 

25 

at  any  but  the  mid-plane. 

For  each  of  the  test  plates,  maximum  advantage  has  been 
taken  of  geometric  symmetries  in  order  to  reduce  the  size  of 
the  finite  element  model.   For  the  tee-stiffened  plate  this 
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presents  no  further  complications  since  the  resulting 
quarter  model  is  non-symmetric.   For  the  plate  with  the 
longitudinally  oriented  rectangular  stiffener,  however,  the 
result  of  half  plate  modelling  is  that  the  stiffener  itself 
lies  on  a  plane  of  symmetry.   This  presents  a  problem  for 
STAGS.   Although  deformations  along  the  symmetric  plane  will 
be  calculated,  no  lateral  displacements  will  be  recognized 
since,  as  far  as  STAGS  is  concerned,  forces  on  either  side 
are  equal.   This  both  affects  the  validity  of  the  displace- 
ments that  are  calculated  and  causes  a  failure  to  indicate 
stiffener  tripping.   STAGS  solves  this  problem  by  allowing 
the  user  to  either  input  initial  geometric  imperfections  in 
the  shell  reference  surface  or  specify  that  the  code 
generate  random  imperfections.   As  stated  in  the  previous 
section,  this  is  done  by  flagging  the  IWIMP  and  I RAMP  vari- 
ables on  the  M-5  record  and  then  following  the  flow  direc- 
tions listed  in  the  manual  through  the  M-7A  record.   For  the 

rectangularly  stiffened  plate  used  here,  imperfections  in 

26 
the  first  three  harmonics  of  0.001  inch  were  specified. 

Finally,  although  not  strictly  related  to  plate  discre- 
tization, there  is  one  other  point  not  mentioned  in  the 
STAGS  user's  manual,  but  which  is  essential  if  the  computer 
code  is  to  successfully  compute  plastic  deformations.   If 


26 

See  the  STAGS1  inputs  for  the  rectangular  stiffened 

plate  in  Appendix  C. 
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plasticity  is  to  be  included,  the  variable  IPLST  on  the  1-1 

record  must  be  set  to  some  value  between  0  and  4  to  indicate 

27 

which  of  four  plastic  strain  theories  should  be  used.     A 

non-zero  positive  value  must  also  be  specified  for  the  vari- 
able NESP.   This  NESP  value  specifies  the  number  of  data 
points  on  a  stress-strain  curve  for  the  plate  material  which 
will  be  input  to  the  1-3  record.   There  are  two  points  that 
the  user's  manual  fails  to  make  clear  about  these  stress- 
strain  points.   First,  this  curve  requires  stress  and  engi- 
neering strains  beyond  the  yield  point.   Second,  the  first 

point  of  the  curve  must  correspond  closely  to  the  modulus  of 

28 
elasticity  of  the  material  specified  in  the  1-2  record. 

All  the  material  and  stress-strain  characteristics  for  the 

6061-T6  aluminum  used  to  manufacture  the  test  plates  here 

were  taken  from  the  Military  Standardization  Handbook, 

Metallic  Materials  and  Elements  for  Aerospace  Vehicle 

Structures  [Ref.  14:   p.  3-184]. 


27 

The  default  value,  0,  automatically  triggers  use  of 

the  Whit e-Besseling  theory,  as  does  the  set  value  1. 

2  8 

That  is,  the  slope  from  the  origin  to  the  first  point 

specified  must  be  the  modulus  of  elasticity.   The  White- 

Besseling  theory  will  apply  approximations  to  attempt  a 

patch-up.   If  the  discrepancy  is  large,  however,  continuity 

between  the  elastic  and  plastic  responses  will  be  lost  and 

the  resulting  output  will  be  in  error. 
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V.   TESTING  THE  CODE 

A.   MODAL  TESTING 

Determination  of  a  structure's  mode  shapes  and  natural 
frequencies  can  provide  a  first  test  of  a  finite  element 
code's  ability  to  predict  the  dynamic  responses  of  an 
excited  system.   Perhaps  of  even  more  importance  to  the 
purposes  of  this  study,  it  can  be  an  excellent  check  of  the 
structural  data  input  to  the  code.   Failure  to  provide 
proper  connectivity  between  shell  units  or  inaccuracies  in 
the  input  dimensions  or  material  properties  will  show  as 
gross  discrepancies  in  the  modal  analysis. 

The  Naval  Postgraduate  School  is  fortunate  to  be  one  of 
the  few  institutions  to  have  a  Hewlett-Packard  System  5451C 
Fourier  Analyzer  with  modal  testing  capability. 

The  technique  used  by  the  H-P  5451C  is  in  many  ways 
similar  to  the  solution  processes  used  by  the  finite  element 
codes  themselves  for  it  involves  solving  basically  the 
same  governing  equation  of  motion  [Ref.  15]. 

Mx  +  Cv  +  K;<  +  f  =  0  (eqn  5.1) 

Since  both  x  an^  f  are  functions  of  time,  the  transfer 
function,  h(x),  may  be  found  by  taking  the  Laplace  trans- 
form of  Equation  5.1 
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B(s)X(s)  =  F(s)  (eqn  5.2) 
where , 

B(s)  =  Ms2  +  Cs  +  K  (eqn  5.3) 
and  , 

H(s)  =  B(s)"1  (eqn  5.4) 

If  the  exciting  forces  and  displacements  are  known,  H(s) 

may  be  readily  found  by  solving  the  equation  which  results 

from  substituting  Equation  5.4  into  Equation  5.2,  namely 

H(s)F(s)  =  X(s)  (eqn  5.5) 

For  an  oscillating  system  of  order  n,  the  transfer  func- 
tion will  have  2n  poles  occurring  in  complex  conjugate  pairs 
Each  pair  of  poles  will  cause  a  mode  of  vibration  in  the 
structure.   Each  of  these  poles  is  a  complex  number  which 
may  be  expressed  as 


Pk  =  -ak  +oui  k  (eqn  5.6) 

where  a,  is  the  damping  coefficient  and  oj,  the  natural 

frequency  for  that  mode. 

Once  a.     and  go,  have  been  determined,  the  resonant 
k       k 

frequency  ft,  (in  radians  per  second)  and  the  damping  ratio 
€,    follow  directly. 
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( eqn  5.7) 

( eqn  5.8) 


The  exciting  forces  and  displacements  that  must  be 
provided  to  the  system  may  be  found  in  a  number  of  ways. 

The  method  chosen  here  was  to  excite  the  plate  using  a 

29  30 

modally  tuned  impulse  hammer. "    An  acceleromet er    was 

attached  to  the  plate  to  record  the  resulting  accelerations, 
from  which  displacements  were  found  by  numerical 
integrat  ion . 

The  two  fundamental  assumptions  of  this  test  procedure 
[Ref .  15 ]  are  that : 

1.  Modal  frequency  and  damping  are  constants  for  any 
transfer  function  taken  from  the  structure. 

2.  Modes  of  vibration  can  be  excited  from  anywhere  on 

an  elastic  structure  except  at  their  node  lines  (where 
no  excitation  is  possible). 

It  should  therefore  be  possible  to  describe  all  modes  of 
interest  by  simply  choosing  a  sufficient  number  of  excita- 
tion points  on  the  structure  and  by  placing  these  points 
close  enough  together  to  minimize  the  possibility  of 
replicat  ion . 

For  the  two  plates  used  in  these  tests  (the  rectangular- 
and  tee-stiffened  test  models),  nine  points  of  excitation 
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30 


PCB  model  086B03. 


PCB  model  302A07 
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31 
were  chosen  on  the  back  of  each  plate.     Care  was  naken  to 

32 

avoid  symmetric  placements  of  the  points.     The  accelerome ler 

was  fixed  at  a  single  location  arbitrarily  chosen  but  spaced 
away  from  the  clamped  boundaries. 

When  any  location  was  tapped  with  the  modal  hammer,  the 
resulting  forces  were  fed  back  through  the  hammer  to  the 
Fourier  Analyzer;  as  were  the  accelerations  from  the  accel- 
erometer.   The  Fourier  Analyzer  then  computed  the  displace- 
ments and  solved  the  Laplace  transform  to  create  the  plot 
shown  in  Figure  5.1.   Each  of  the  peaks  in  this  representa- 
tion corresponds  to  the  poles  of  the  transfer  function  for  a 
single  point  of  excitation.   By  then  combining  the  plots  for 
all  points,  an  accurate  representation  of  the  system  modes 
was  obtained.   It  was  then  simply  a  matter  of  analyzing  each 
of  the  poles  in  descending  order  of  magnitude  to  whatever 
total  number  of  modes  is  of  interest  to  the  user  to  deter- 
mine the  dominant  modes  of  the  system. 

B.   UNDERWATER  SHOCK  TESTING 

Underwater  explosion  testing  for  all  studies  conducted 
in  this  series  to  date  has  been  done  at  the  West  Coast  Shock 
Test  Facility  (WCSF),  located  at  the  old  Hunter's  Point 
Naval  Shipyard  near  San  Francisco,  California.   This  is  a 


31 

The  plates  were  placed  face  downward  upon  a  solid 

marble  slab. 

32 

The  actual  number  and  spacing  of  the  excitation  points 

seems  to  be  primarily  a  matter  of  experience. 
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Figure  5.1   Typical  Modal  Test  Plot 
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Navy  activity  reporting  to  the  Supervisor  for  Shipbuilding, 
Conversion,  and  Repair,  San  Francisco.   As  the  only  licensed 
facility  of  its  type  on  the  West  Coast,  they  have  gained  a 
considerable  amount  of  experience  in  conducting  underwater 
explosion  testing;  although  most  of  it  has  been  large  scale 
work  in  support  of  the  Navy's  shock  qualification  program 
rather  than  experimentation. 
1 .   The  Test  Platform 

In  order  to  provide  the  air  backed  condition  that 
was  required  to  simulate  a  ship's  hull  and  to  ensure  rigidly 
clamped  boundary  conditions,  the  test  plates  used  in  this 
series  of  studies  were  securely  bolted  to  a  heavy  steel 
backing  structure  (Fig.  5.2).   This  structure  consisted  of  a 
box  constructed  of  lj    inch  thick  A6  structural  steel.   A 
wide  flange  was  welded  around  the  open  side  of  the  box  to 
which  the  edge  of  the  plate  was  bolted.   An  O-ring  .gasket 
was  fit  into  a  channel  machined  into  the  surface  of  the 
flange  to  provide  watertight  integrity.   Access  to  instru- 
mentation on  the  inside  of  the  box  was  allowed  by  a  subma- 

33 
rine  style  penetrator  affixed  through  the  bottom. 

To  perform  the  actual  testing,  the  plate  and  backing 

structure  were  suspended  by  heavy  steel  cables  below  two 

large  pneumatic  floats  in  the  configuration  developed  in 


33 

For  a  more  complete  description  of  this  backing  struc- 
ture and  the  rationale  used  in  designing  it,  see  Reference 
3,  pp.  88-94. 

90 


~'f.  ■>i?y^S?r~. 


Figure  5.2   Test  Plate  and  Backing  Structure 
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Chapter  2  (Fig.  5.3).   The  cables  were  cut  to  length  to 
assure  that  the  nearest  edges  of  the  test  panel  stiffeners 
were  maintained  at  depths  as  close  to  thirteen  feet  below 
the  surface  as  possible.   Once  the  entire  apparatus  had  been 
towed  into  the  bay,  the  explosive  charge  was  lowered  into 
position  nine  feet  above  the  stiffeners  in  line  vertically 
with  the  center  of  the  plate. 
2 .   Instrumentation 

In  his  work,  which  has  already  been  cited  several 
times,  LT  Rentz  outlined  the  basic  theory  governing  the  test 
of  the  flat  plate  and  began  developing  the  procedures  which 
will  be  enlarged  upon  here. 

Essentially,  the  information  sought  for  comparison 
with  the  computer  simulation  codes  being  tested  has  been: 

1.  Fluid  pressures  at  the  plate  and  in  the  freefield. 

2.  Strain  in  the  plate. 

3.  Final  deflections  in  the  plate. 
Instrumentation  used  to  gather  this  data  was  as  follows: 

a.   Pressure  Measurements 

Three  pressure  gages  were  used  in  the  tests 
conducted  here.   Two  of  these  were  located  just  as  they  had 
been  in  the  previous  study  [Ref.  3:   p.  46].   One  was  located 
one  foot  above  the  center  of  the  plate  and  the  other  was 
clamped  to  the  aluminum  block  bolted  to  the  side  of  the 
plate.   The  forward  edge  of  this  second  pressure  gage  was 
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Figure  5.3   Test  Platform 
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positioned  so  as  to  be  on  an  arc  of  the  same  radius  as  the 
standoff  distance. 

The  pressure  recordings  that  had  been  obtained 
from  these  two  gages  had  been  so  subject  to  interference 
from  the  plate  structure  and  the  rigging  of  the  test  plat- 
form that  the  futility 'of  conducting  further  tests  without  a 
true  free  field  pressure  gage  was  readily  obvious.   Accor- 
dingly, although  the  first  two  gages  were  retained  as 
backups,  a  third  gage  was  positioned  on  an  arc  equal  to  the 
standoff  distance  but  removed  five  feet  from  the  centerline 
of  the  plate  (Fig.  5.4).   As  will  be  discussed  in  the 
following  chapter,  the  results  obtained  from  this  gage 
during  the  test  shot  on  the  tee-stiffened  plate  were 
excellent . 

All  tests  in  this  series  utilized  one-quarter 

inch  tourmaline  crystal  pressure  gages  rated  to  10,000  psi. 

34 

The  response  ratio,    R  ,  is  a  function  of  the  decay 

constant  for  the  charge,  9,  divided  by  the  transit  time  of 
the  shock  wave  across  the  gage,  tQ  [Ref.  16].   For  the 
eight  pound  TNT  charge  chosen  as  the  standard  for  these 
studies , 

9  =  0.15  msec 


34 

Defined  as  the  ratio  of  the  apparent  peak  pressure  to 

the  actual  peak  pressure. 
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Figure    5.4      Location    of    Free    Field    Pressure    Gage 


95 


For  a  quarter  inch  gage, 
0.25  in 


t 


60,00c1 


sec 


=  4.167  x  10  6  sec 


Applying  the  correction  factor  for  the  oil-filled  protective 
boot  which  surrounds  the  crystal, 

1.5  tD  =  6.25  x  10~6  sec 

Thus  , 


=  24 


1.5  tD 


From  the  chart  on  page  17  of  Reference  16. 

R   =  0.9733 

P 

This  indicates  that  a  very  good  correlation  can  be  expected 
between  the  actual  and  measured  pressures  ,  and  is  also  the 
best  indication  that  the  correct  gage  size  is  being  used, 
b.   Strain  Measurements 

Perhaps  the  single  most  critical  aspect  of 
obtaining  successful  data  from  underwater  explosion  testing 
of  the  type  conducted  here  is  the  choice  and  application  of 
the  strain  gages. 
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For  the  tee-stiffened  plate,  twenty  SR-4  type 

35 
FAE-25-12-S13ET  strain  gages   "  were  attached  as  shown  in 

Figure  5.5.   These  were  the  same  type  used  in  Reference  3; 
and,  although  rated  to  only  2000  microstrain,  it  had 
been  demonstrated  that  their  range  could  be  extended  suffi- 
ciently to  provide  useable  data  over  the  very  earliest 
increments  of  the  shock  response.   As  in  the  earlier  test, 
BLH ' s  SR-4  Permabond  910  was  used  as  the  fixative,  and  the 
attached  strain  gages  and  their  wiring  was  covered  with 

PR-1422  A-l  Base  Compound,  manufactured  by  Products  Research 

36 
and  Chemical  Corporation  of  Glouster  City,  NJ .     This  base 

compound  was  also  used  on  the  plate  with  the  rectangular 

O  '" 7 

stiffener  where  19  MX  type  PAHE-250BA-350EN  strain  gages 
(Fig.  5.6)  rated  to  approximately  200,000  microstrain  were 
used  in  an  attempt  to  increase  the  time  span  over  which 
information  could  be  gathered.   Also,  at  the  recommendation 

of  the  distributor,  BLH  type  EPY-150  epoxy  adhesive  was 

38 
used  instead  of  the  Bondo  910. 
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Manufactured  by  BLH  Electronics  of  Walt  ham,  MA 


36 

This  was  done  to  provide  and  smooth  compliant 

covering  that  would  keep  the  gages  and  wiring  from  being 
torn  off  with  the  passage  of  the  shock  wave,  and  to  seal  the 
gages  from  saltwater  once  the  rig  was  submerged. 


37 


Manufactured  by  Experimental  Stress  Analysis  Products. 


O  Q 

In  theory,  the  epoxy  provides  a  stronger  bond  than 
does  the  Bondo,  which  is  a  contact-type  adhesive. 
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In  practice,  although  some  useful  data  was 
obtained  from  each  test,  neither  was  a  particularly 
successful  source  of  strain  informal  ion .   The  number  of 
gages  lost  far  exceeded  the  number  that  performed. 

In  recapitulation,  a  number  of  faults  were 
found  for  which  the  following  corrective  measures  have  been 
proposed . 

1.  Although  the  MX  type  strain  gages  greatly  exceed 
the  strain  requirements  expected  in  these  tests, 
their  use  is  strongly  recommended.   Little  or  no 
loss  of  accuracy  is  experienced  in  the  lower  strain 
regions  through  use  of  the  high-elongation  gages, 
whereas  mid-range  gages  tend  to  peak  right  at  the 
strain  ranges  expected. 

2.  On  any  strain  gage,  the  weakest  points  are  the  sold- 
ering tabs  to  which  the  gage  wires  are  attached. 
These  are  brittle  and  will  break  if  any  bending 
stress  is  applied.   All  motion  should  be  restricted 
to  the  plane  of  the  gage.   This  means  that  the  strain 
gage  wires  must  be  extremely  flexible  so  that  any 
bending  motion,  either  during  installation  or  as  a 
result  of  plate  deflections,  is  not  transmitted  to 

the  gage  tabs.   Rated  strain  gage  wire*  is 

39 

recommended . 


39 

On  the  tee-stiffened  plate,  28  gage  wire  made  specifi- 
cally for  strain  gages  was  used.   This  proved  to  be  quite 
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3.   Even  with  flexible  strain  gage  wire,  a  determined 

force  along  the  length  of  the  wire,  as  would  be  expert 
enced  during  plate  deformation  will  break  the  tab, 
the  wire,  or  the  solder  connection  between  the  two. 
Although  a  plastic  medium,  the  PRC  Base  Compound 
used  to  protect  the  gages  and  wires  was  not  flexible 
enough  to  allow  the  wires  to  expand  with  the  plate. 
A  suggested  solution  to  this  problem  is  to  lay  the 
gage  wires  with  some  s-curves  to  allow  elongation  as 
the  plate  deforms.   The  wires  should  not,  of  course, 
be  otherwise  affixed.   A  piping  of  Barrier  WD,  a 
tacky,  semi-fluid  substance  manufactured  by  BLH  which 
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looks  somewhat  like  a  combinat ion  of  beeswax  and 
petroleum  jelly,  could  then  be  applied  over  both  gages 
and  wires.     This  piping  of  Barrier  WD  should  then  be 
covered  by  a  layer  of  aluminum  foil  and  then  the  whole 
overlayed  by  the  usual  PBC  Base  Compound.   The  net 
effect  is  an  armored  "tube"  of  Barrier  WD  through  which 
the  wires  can  run  when  the  plate  deforms. 
The  fixative  used  to  attach  the  strain  gages  should, 
itself,  be  a  subject  for  careful  consideration.   Even 
though  the  epoxy  is  supposed  to  be  the  stronger  medium, 
experience  gained  here  indicates  that  it  is  not  as 
satisfactory  as  the  Bondo  for  explosion  testing.   The 
postulated  reason  for  this  seeming  discrepancy  is  that 
the  epoxy  places  a  thicker  layer  of  adhesive  between 
gage  and  plate  material.   Given  the  extremely  rapid 
loadings  experienced,  intense  shears  are  created  in  the 
thicker  substance  and  the  epoxy  literally  cracks  and 

spalls  off.   The  Bondo  910  therefore  seems  to  be  the 

42 
preferable  adhesive. 


41 

A  bakers  pastry  tube  works  very  nicely  to  do  this. 

42 

No  failures  due  to  loading  were  noted  in  the  strain 

gages  attached  with  Bondo.   Therefore,  Even  if  the  cracking 
problem  had  not  occurred  in  the  epoxy,  the  greater  separation 
between  gage  and  plate  material  that  the  epoxy  creates  would 
argue  in  favor  of  the  Bondo. 
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c.   Plate  Deflections 

Plate  deflections  were  by  far  the  most  reliable 
source  of  data  that  could  be  obtained  as  a  result  of  testing 
since  no  electronics  were  involved.   To  get  these,  precise 
measurements  were  made  at  specified  points  on  the  plate  both 
before  and  after  the  test  shot.   A  jig  was  devised  which 
took  advantage  of  the  rigid  plate  boundaries  to  provide  a 
reference  table  for  a  dial  indicator  depth  gage  (Fig.  5.7). 
All  measurements  were  made  while  the  test  plate  was  bolted 

onto  the  backing  structure  to  ensure  that  plate  warpage  did 

43 

not  affect  the  readings. 

C.   DATE  RETRIEVAL 

Signals  from  the  gages  mounted  on  and  about  the  test 
plate  were  fed  back  through  shielded  cable  to  a  common 
floating  terminal  box.   From  there,  they  were  passed  to  an 
instrumentation  shack  located  on  one  of  the  jetties  defining 
the  small  bay  in  which  the  testing  was  conducted.   A  sche- 
matic of  the  setup  is  shown  in  Figure-  5.8.   Calibration 
levels  were  applied  to  all  pressure  signals  by  Endevco 
Model  4470  Signal  Conditioners.   These  same  signal  condi- 
tioners were  used  to  provide  calibration  for  all  but  twelve 


\ 
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Deflections  have  already  been  given  for  the  rectan- 
gular stiffened  plate  in  Figures  3.7  and  3.8.   No  post-shot 
measurements  were  made  of  the  tee-stiffened  plate  because  of 
the  heavy  damage  which  was  inflicted  when  the  plate  was  hit 
by  one  of  the  supporting  cables. 


103 


Figure  5.7   Jig  Used  to  Measure  Plate  Deflections 


104 


PRESSURE 
GAGES 

V 


FLOATING 
TERMINAL  BOX 


*C 


STRAIN 
GAGES 


LANDLINE  CABLES 


VOICE 


INSTRUMENT  ROOM 

-TERMINAL  BOARD 


r 


i 

JL 


FOURIER 
ANALYZER 


FIRST  12  CHANNELS 


BRIDGE 
BALANCE 


DC 
POWER 


CHANNELS 
7-13-21 


SIGNAL 
CONDITIONER 


AMPLIFIERS 


f 


FM  TAPE 

RECORDER 


I 


TIME 


1  CODE 


i 
PLAYBACK 

T 


AMPLIFIERS 


VISICORDER 


OSCILLO- 
SCOPE 


Figure  5 . 8 


Flow  Chart  for  Electronics  Set-Up  at  WCSF 


105 


of  the  strain  gage  channels.   Those  remaining  twelve 
channels  were  calibrated  by  two  six-channel  Honeywell  Model 
B2-6  Universal  Bridge  Balance  units.   The  strain  gage  signals 
were  then  amplified  by  Bell  and  Howell  CEC  1-168  DC  ampli- 
fiers before  final  recording.   All  signals  were  recorded 
on  two  twelve  channel  Amp.ex  FR  1300  tape  recorders. 

Calibration  levels  for  the  tourmaline  crystal  pressure 
gages  were  all  set  at  the  10,000  psi  theoretical  maximum. 
For  the  strain  gages,  however,  attempts  were  made  to  use  the 
computer  codes  to  provide  estimates  of  expected  strains 
using  the  pressure  and  decay  constant  calculated  from 
Equations  2.2  and  2.3.  This  was  especially  critical  for 
the  strain  levels  for  the  tee-stiffened  plate  where  some  of 
the  strain  gages  were  expected  to  be  operating  at  or  beyond 
their  rated  limits. 

Unfortunately,  difficulties  encountered  in  learning 

to  use  the  USA-STAGS  code  prevented  runs  from  being -made 

44 
before  the  shot  dates. 

Advantage  could  have  been  taken  of  earlier  estimates 
used  for  the  plate  tested  in  Reference  3  if  the  charge 
weights  had  remained  the  same;  however,  such  was  not  the 
case  . 

Since  8  pound  TNT  charges  are  listed  as  Navy  standard 
stock  items  and  order  had  been  placed  well  in  advance,  no 
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In  fact  no  successful  runs  were  made  until  long  after 


& 


all  experimentation  had  been  completed 
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problems  were  anticipated  in  obtaining  charges  of  this  size. 
As  the  date  of  the  first  test  neared ,  however,  it  became 
obvious  that  no  8  pound  charges  would  be  available.   The 
nearest  alternate  that  could  be  provided  in  time  for  the 
scheduled  test  was  a  9i    pound  TNT  charge. 

It  can  perhaps  be  appreciated  that  dimensional  changes 
to  the  test  platform  are  not  simply  a  matter  of  moving  the 
charge  location.   As  charge  weights  are  changed,  the  stand- 
off distance  and  the  charge  depth  both  must  be  changed; 

which  involves  cutting  new  supporting  cables  and  recali- 
brating flotation  depths.   All  of  which  is  extremely  time 
consuming.   When  the  personnel  at  WCSF  were  contacted  about 
the  possibility  of  changing  the  configuration  for  a  9f 
pound  charge,  they  advised  that  it  could  not  be  done  in  time 
for  the  scheduled  tests. 

The  choice  for  the  first  test  was  clear,  either  run  the 
test  with  the  existing  platform  and  a  9?  pound  charge,  or 
abort  the  experiment .   It  was  determined  that  more  benefit 
could  be  gained  from  conducting  the  test. 

The  values  for  Pn  and  6  found  from  Equations  2.2  and 
2.3  for  a  9-|  pound  charge  are: 


P   =  4081.7  psi,  9  =  0.161  msec 


This  initial  pressure  was  markedly  higher  than  that  for 
an  8  pound  charge.   The  maximum  bubble  radius,  too,  was 
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much  larger.   From  Equation  2.6,  A     =  14.03  ft.   This 

°  max 

gives  a  charge  depth  to  bubble  radius  ratio  of 


4  ft 

=  0.285 


14.03  ft 


which  implies  that  some  pressure  release  will  take  place 
into  the  atmostphere. 

For  the  second  test  shot  ,  the  problem  was  quite 
different.   As  was  discussed  in  Chapter  III,  the  tee- 
stiffened  model  appeared  too  rigid  to  successfully  demon- 
strate stiffener  tripping.   The  plate  with  the  longitudinal 
rectangular  stiffener  was  a  largely  tentative  attempt  to 
produce  that  action.   It  was  desired  to  view  the  responses 
of  the  stiffener  under  a  variety  of  loadings,  beginning  with 
small  charges  and  concluding  with  ones  in  the  8  to  9|  pound 
range.   Again  use  was  made  of  what  was  available.   In  this  " 
case,  the  available  charge  weight  was  h    pound. 

For  a  \    pound  charge,  the  maximum  bubble  radius  is  3.02 
feet;  which  leads  to  a  charge  depth  to  bubble  radius  ratio 
of  1.33.   This  implies  that  bubble  pulsation  will  be  a 
problem.   Fortunately,  however,  the  time  of  the  pulse  is 
relatively  late,  t  =  0.17  seconds.   Thus,  by  limiting  the 
period  of  study  to  early  times  the  bubble  pulse  can  be 
effectively  neglected. 

Calibration  levels  for  the  test  shot  on  the  tee- 
stiffened  plate  were  determined  by  making  an  EPSA  run  on  the 
plate  model  used  in  Reference  3  but  with  a  9|  pound  charge 
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weight.  As  has  already  been  stared,  geometric  similarities 
between  these  two  plates  permitted  broad  generalizations  to 
be  made  between  them  even  if  exact  comparisons  could  not  be 
drawn.   The  results  of  that  run  are  shown  in  Table  III. 

No  similar  comparisons  could  be  made  with  the  plate 
with  the  longitudinal  rectangular  stiffener.   Strain  gage 
calibrations  were  therefore  arbitrarily  set  at  a  value  at 
the  higher  end  of  the  strain  gage  range  which  could  be 
expected  to  exceed  all  requirements.   In  this  case,  that 
value  was  150,000  microstrain. 

The  results  of  these  tests  will  be  discussed  in  the 
next  chapter. 

Table  III 
Calibration  Data  for  Tee-Stiffened  Plate 

Gage  EPSA  Strain*  Gage  EPSA  Strain* 

1  6210  11 

2  6210  12 

3  17600  13 

4  8500  14 

5  21700  15 

6  69400  16 

7  69400  17 

8  69400  18 

9  100000  19 
10               49500              20 

*A11  values  in  microstrain. 


76800 

7610 

7610 

8510 

7610 

21700 

6210 

6210 

17600 

58200 
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VI.   RESULTS  AND  CONCLUSIONS 

A.   MODAL  COMPARISONS 

As  stated  in  chapter  four,  the  primary  purpose  of  modal 
testing  is  to  provide  a  check  of  the  finite  element  model. 
A  properly  constructed  model  should  exhibit  substantially 
the  same  fundamental  frequency  as  that  determined  by  modal 
testing  and  should  demonstrate  close  coincidence  with  other 
frequency  values  across  the  range  tested.   Exact  concurrence 
cannot  be  expected  since  the  computer  code  is  analyzing  an 
idealized  model  which  will  not  be  affected  by  the  externali- 
ties inherent  in  the  actual  test.   Also,  the  number  of 
frequencies  determined  by  the  computer  is  a  function  of  the 
number  of  iterations  than  can  be  made  over  a  given  time 
span.   Certain  economies  in  the  use  of  computer  time  must 
therefore  be  observed.   While  the  code  might  be  able  to 
calculate  the  frequencies  for  a  large  number  of  modes;  such 
calculation  will  be  purchased  with  a  heavy  expenditure  in 
CPU  time  which  is  not  really  necessary  for  an  adequate 
check . 

Accordingly,  the  scheme  followed  here  has  been  to  limit 
the  number  of  modal  frequencies  computed  by  the  code  to 
those  that  could  be  determined  in  two  runs  for  each  plate. 
Each  of  these  runs  can  be  performed  entirely  within  the 
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45 
combined  STAGS1 /STAGS2    processor  without  reference  to  the 

USA  add-on.   The  first  run  simply  specifies  a  broad  range 

within  which  the  computer  can  search  for  eigenvalues;  the 

second  uses  one  of  the  mid-range  frequencies  determined  from 

the  first  run  eigenvalues  as  a  "seed"  for  further  iterations 

46 
in  order  to  produce  a  refined  product  .     The  actual  perform- 
ance of  these  functions  is  controlled  by  the  D-2  and  D-3 
records  ofthe  STAGS  input  deck.   Since  the  procedure  is 
outlined  clearly  in  Reference  13  and  example  input  files  are 
provided  for  both  the  tee-stiffened  and  rectangular  stiff- 
ened plates  in  Appendix  C,  further  elaboration  will  not  be 
provided  here. 

Modal  testing  of  the  milled  plates  was  conducted  in  the 
manner  described  in  Section  D  of  Chapter  IV.   The  results  of 
both  that  testing  and  the  STAGS2  estimates  are  provided  for 
comparison  in  Tables  IV  and  V  tee-stiffened  and  rectangular- 
stiffened  plates,  respectively. 

It  will  be  noted  that,  for  each  of  these  plates,  the 
STAGS2  predicted  values  seem  to  correspond  very  well  with 


45 

Generally  referred  to  here  as  a  STAGS2  run. 

46 

Since,  in  this  case,  modal  testing  of  the  milled 

plates  was  conducted  prior  to  the  computer  runs,  essentially 
the  same  frequency  range  was  used  for  the  computer  runs  as 
for  the  testing.   Over  this  range,  the  first  STAGS2  run  will 
compute  five  modes  with  their  corresponding  natural  frequen- 
cies.  Although  the  user's  manual  states  that  calculation  of 
up  to  twenty  eigenvalues  can  be  specified,  and  that  all 
existing  eigenvalues  up  to  that  number  will  be  computed, 
five  seems  to  be  the  true  maximum.   No  increase  in  the 
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Table  IV 


Modal  Comparisons  for  Tee-Stiffened  Plate 


Measured  Resonant 
Frequency  (Hz-) 

629.7410 
653.2053 
721.1747 
772.3517 
.  793.0253 
1022  .4708 
1310.2979 
1447.0071 
1448.4597 
1480.9890 
1499.5906 
1551.0845 
1556.0732 
1590.2947 
1747.9771 
1994.4626 
2143.0918 
2236.9248 
2259.6445 
2275.5000 
2430.0249 
2472 .2305 
2571.7832 


STAGS2  Predicted 
Modal  Frequency  (Hz) 


772.7611 


957.8445 


2043.4048 
2200.9004 


2629.4729 
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Table  V 
Modal  Comparisons  for  Rectangular-Stiffened  Plate 


Measured  Resonant 
Frequency  (Hz) 


516 

5652 

804 

2715 

924 

.8018 

980 

5009 

981 

1315 

1001 

6130 

1118 

0952 

1204 

4885 

1229 

9939 

1365 

3230 

1428 

2920 

1651 

.2498 

1651 

.2747 

1852 

.1162 

2040 

.0291 

2154 

.0908 

2368 

.5132 

2603 

.9736 

2725 

.8828 

STAGS2  Predicted 
Modal  Frequecy  (Hz) 


790.0587 


1115.3547 


1812 .0698 
2048.0015 
2176 .5366 


3039.6057 
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values  from  the  modal  testing,  the  greatest  relative  error 
(using  the  modal  test  value  as  a  base)  being  6. 32"  for  the 
tee-stiffened  plate  and  2.16%  for  the  rectangular-stiffened 
plate,  and  the  average  relative  errors  being  3.16%  and  1.12% 
for  the  tee-  and  rectangular-stiffened  plates,  respectively. 

Less  promising  is  the  seeming  failure  of  the  code  to 
predict  the  fundamental  frequency  for  either  of  the  plates. 
It  is  suggested,  however,  that  this  is  not  due  to  any  errors 
in  the  code  or  in  the  inputs  to  the  code  ,  but  rather  is 
attributable  to  the  fact  that  the  boundary  conditions 
actually  tested  were  not  the  clamped-clamped  conditions 
specified  in  the  STAGS  inputs.   The  H-P  5451C  Fourier 
Analyzer  appears  to  be  quite  sensitive  to  relatively 
small  changes  in  the  boundary  conditions.   When  the  tests 
were  conducted,  the  plates  were  laid  face-downward  upon  a 
heavy  marble  slab.   No  attempt  was  made  to  further  clamp 
the  heavy  aluminum  frame  into  which  the  test  plates  are 


amount  of  CPU  time  allowed  for  eigenvalue  analysis  (NSEC  on 
the  D-2  record)  seems  to  affect  that  number.   Computation  of 
eigenvalues  over  a  range  without  any  seeding  of  known  values 
is  a  slow  process  prone  to  convergence  problems.   It  should 
therefore  only  be  used  to  derive  a  value  which  can  be  used 
for  the  second  run;  or  it  can  be  eliminated  altogether  if 
one  of  the  natural  frequencies  is  already  known.   If  such  a 
frequency  has  been  determined,  it  can  be  entered  directly 
into  the  code  and  all  frequencies  closest  to  multiples  of  a 
fixed  interval  (SHIFT  on  the  D-3  record)  will  be  computed. 
This  process  is  reasonably  fast  and  will  produce  all 
existing  frequencies  up  to  the  number  (NEIG)  specified. 
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milled,  nor  was  it  possible  to  ensure  an  absolutely  flush 
mount  between  plate  and  slab. 

When  the  tee-stiffened  plate  was  rested,  the  bolt  holes 

had  already  been  drilled  in  the  frame.   These  bol"C  holes 
seemed  to  cause  reflections  which  were  picked-up  by  the 
Fourier  Analyzer  as  additional  modes.   These  may  be  observed 
to  some  degree  in  Table  IV  in  the  repeated  tendency  for  a 
group  of  frequencies  to  cluster  about  a  single  value.   The 
effect  was  far  more  noticeable  on  the  oscilloscope  traces 

obtained  as  each  of  the  excitation  points  on  the  plate  was 

47 
loaded.     When  all  of  the  variables  are  considered,  it  is 

believed  that  the  actual  correspondence  between  the  modal 

test  results  and  the  STAGS2  prediction  is  quite  good  and 

provides  adequat-e  confirmation  of  both  the  code's  ability 

to  predict  the  natural  frequencies  of  a  simple  structure 

and  of  the  finite  element  models  input  to  the  code  for  this 

study  . 

B.   COMPARISONS  WITH  UNDERWATER  TEST  RESULTS 

As  has  already  been  stated  in  Chapter  V,  all  underwater 
shock  testing  for  this  series  of  studies  was  conducted  at 
the  West  Coast  Shock  Test  Facility.   The  outputs  from  the 


47 

It  is  strongly  recommended  that  future  modal  testing 

be  conducted  before  the  bolt  holes  are  drilled  in  the  frame 
An  attempt  might  also  be  made  to  weight  the  frame  with  sand- 
bags to  simulate  as  nearly  as  possible  the  clamped  boundary 
conditions  input  to  the  code. 
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pressure  and  strain  gages  used  for  the  tests  along  with 
their  associated  calibration  and  time  signals  were  recorded 
by  two  Ampex  FR  1300  tape  recorders.   After  each  of  the  test 
shots,  the  data  tapes  and  one  of  the  tape  recorders  were 
transported  to  the  Naval  Postgraduate  School  where  the  tapes 
were  played-back,  filtered,  and  recorded  on  disk  on  the  H-P 
5451C  Fourier  Analyzer  for  later  analysis  and  plotting. 

As  thorough  a  pretest  analysis  as  was  possible  without 
actual  disassembly  was  conducted  upon  the  equipment  at  the 
West  Coast  Shock  Test  Facility  before  each  shot.   Calibra- 
tion levels  for  all  sensing  devices  were  also  checked 
using  technicians  and  equipment  provided  by  the  Naval 
Postgraduate  School.   Despite  these  precautions,  the  quality 
of  data  obtained  from  the  tests  was  highly  problematical. 

An  excellent  free-field  pressure  history  was  obtained 
from  the  test  shot  conducted  on  the  tee-stiffened  plate 
(Fig.  6.1);  fully  compensating  the  difficulties  encountered 
in  mounting  the  third  pressure  gage  in  a  position  where  it 
would  not  be  affected  by  interference.   As  can  readily  be 
observed,  this  curve  displays  all  of  the  classical  exponen- 
tial decay  characteristics  predicted  in  Chapter  II.   Since 
this  is  so,  the  decay  constant  will  correspond  approximately 

to  the  time  it  takes  for  the  pressure  to  drop  to  —  (about 

e 

.36788)  of  its  peak  value.   From  Figure  6.1,  the  peak 
incident  pressure  was  3125  psi.   Therefore, 
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Figure  6.1   Freefield  Pressure  for  Tee-Stiffened  Plate 


117 


(3125  psi)(. 36788)  =  1149.6  psi 

which  occurs  at  about  1.56  milliseconds. 

These  values  were  input  to  the  USA-STAGS  code  to  conduct 
a  post-shot  analysis  for  comparison  to  test  data. 

The  peak  pressure  and  decay  constant  obtained  from  the 
test  results  are  much  smaller  than  those  predicted  by 
Equations  2.2  and  2.3.   In  retrospect,  this  is  not  parti- 
cularly surprising.   First,  it  will  be  recalled  that  Equations 
2.2  and  2.3  are  derived  from  empirical  results  and  are  there- 
fore subject  to  wide  variations.   Even  more  pertinent,  it 
was  found  in  Chapter  V  that  the  ratio  of  the  charge  depth  to 
the  maximum  bubble  radius  was  much  less  than  the  .50 
threshold  for  pressure  loss.   The  decreased  pressure  and 
more  gradual  decay  are  exactly  what  would  be  expected  from  a 
charge  that  had  vented  a  significant  portion  of  its  energy 
in  the  first  moments  after  detonation. 

Given  the  good  results  of  the  first  shot,  a  high  degree 
of  optimism  was  established  for  the  results  of  the  pressure 
data  from  the  second,  the  k    pound  charge  on  the  plate  with 
the  longitudinal  rectangular  stiffener.   That  optimism  was 
not  justified. 

Upon  playback  of  the  data  for  the  second  shot,  it  became 
clear  that  the  amplifiers  from  all  three  of  the  pressure 
gages  had  become  totally  saturated.   This  was  surprising 
since  the  pressure  gages  were  calibrated  for  10,000  psi  and 
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the  peak  pressure  for  a  §  pound  TNT  charge  is  far  less 
than  that . 

The  methods  by  which  the  two  tests  were  conducted 
suggested  a  probable  cause  for  the  pressure  gage  saturation. 

Once  all  wiring  was  connected  and  the  gages  tested,  the 
test  platform  was  lowered  into  the  waters  of  the  bay.   Heavy 
halyards  were  attached  to  either  end  of  the  floats,  one  to 
be  used  as  an  outhaul  and  the  other  as  an  inhaul  .   Using 
these,  the  platform  was  positioned  in  the  bay  a  safe 
distance  from  the  jetties  before  the  explosive  charge  was 
installed  and  the  test  conducted. 

To  conduct  the  test  with  the  9k    pound  charge,  the 
platform  was  centered  in  the  bay  to  keep  it  as  far  away  from 
the  jetties  as  possible.   There  was  a  heavy  wind  blowing  and 
when  the  cable  used  to  carry  the  firing  current  to  the 
explosive  charge  was  taken  out  to  the  platform  it  was  blown 
well  away  from  the  platform  and  the  instrument  cables. 
Since  it  was  separated  from  the  other  instruments,  the  free- 
field  pressure  gage  cable  was  farthest  of  all  from  the 
detonating  cable. 

When  the  test  shot  was  conducted  on  the  rectangular- 
stiffened  plate,  the  charge  was  sufficiently  small  that  the 
platform  was  kept  much  nearer  pierside  to  make  it  easier  to 
handle.   All  of  the  cabling  was  run  out  to  the  platform 
along  the  inhaul  . 
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The  voltage  used  to  detonate  the  explosive  charge  is 
quite  large  (800  volts  DC)  and  sets  up  a  noticeable  induc- 
tance on  all  electrical  cabling  in  the  vicinity.   In  fact, 
on  the  strain  gages,  which  are  very  low  impedance  devices 
and  therefore  relatively  insusceptible  to  inductive  noise 

this  was  used  as  a  convenient,  if  not  very  accurate,  zero 

,48 
time  mark. 

The  pressure  gages,  on  the  other  hand,  are  extremely 
high  capacitance,  high  impedance  gages  and  correspondingly 
are  very  much  affected  by  induced  currents.   It  is  therefore 
believed  that  the  close  proximity  of  the  detonation  cable 
was  the  sole  cause  of  the  saturation  observed  on  the  pres- 
sure gages  during  the  second  test.   Be  that  as  it  may,  no 
pressure  history  was  obtained  for  the  rectangular-stiffened 
plate,  making  it  necessary  to  use  the  peak  pressure  and 
decay  constant  calculated  from  Equations  2.2  and  2.3  in 
the  subsequent  USA-STAGS  run. 

The  lack  of  a  zero  time  signal  was  not  the  only  defi- 
ciency noted  on  the  strain  gage  histories.   Three  others 
should  also  be  mentioned  by  way  of  warning  for  those 
conducting  similar  studies. 


48 

This  was  a  major  deficiency  in  the  way  testing  was 

conducted  at  WCSF.   No  accurate  zero  time  was  provided  on 

any  channel.   The  detonation  signal  was  used  as  a  best 

alternative,  but  since  that  typically  covered  a  span  of 

about  one  millisecond  it  was  almost  useless  when  only  the 

first  few  milliseconds  were  being  studied. 
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Before  the  test  platform  was  placed  in  the  water,  all  of 
the  strain  gage  signals  were  balanced  with  their  calibration 
resistors,  setting  their  signals  to  zero.   As  soon  as  the 
platform  was  placed  in  the  water,  the  change  in  temperature 
effected  a  corresponding  change  in  the  resistance  of  the 
strain  gages.   This  in  turn  changed  the  balance  of  the 
bridge.   This  was  an  ongoing  problem  throughout  the  period 
of  the  tests;  many  of  the  strain  gages  never  were  brought 
into  proper  balance.   The  problem  was  exacerbated  by  the 
fact  that  the  strain  gage  and  its  associated  wiring  formed 
an  extremely  long  leg  of  the  balance  bridge  which  was 
subject  to  the  environment  of  the  bay.   The  balancing  leg 
consisted  of  a  calibration  resistor  located  in  the  instru- 
ment shack;  a  very  different  environment  from  that  outdoors. 

The  second  of  these  three  problems  is  one  which  is 
readily  correctable  by  the  experimenter  if  he  is  aware  of 
it.   When  the  test  plate  and  backing  structure  were  sent  to 
the  West  Coast  Shock  Test  Facility,  all  strain  gage  leads 
were  connected  to  junctions  in  the  terminal  box  located  on 
the  side  of  the  backing  structure.   Connecting  these  leads 
through  long  cable  runs  was  one  of  the  many  tasks  that  was 
performed  by  WCSF  personnel  between  the  time  that  the  plate 
was  delivered  and  the  date  of  the  test  shot.   Their  experi- 
ence was  primarily  with  the  Navy's  shock  test  program, 
wherein  the  magnitudes  of  the  strain  outputs  were  of  more 
value  than  the  sign.   Therefore,  little  attention  was  paid 
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to  the  polarity  of  the  strain  gage  connections.   Thus,  the 
histories  recorded  can  be  the  exact  inverse  of  the  true  his- 
tory.  This  can  be  especially  debilitating  when  comparisons 
of  the  shapes  of  strain  curves  is  of  particular  importance. 

Even  more  debilitating,  but  far  less  obvious,  was  the 
last  of  the  problems  noted. 

The  amplifiers  on  the  Ampex  FR  1300  tape  recorders  used 
to  record  the  strain  signals  are  calibrated  for  a  40°o  devia- 
tion from  a  one  volt  rms  input  signal,  peak-to-peak.   The 
maximum  signal  deflection  that  is  theoretically  possible  is 
about  1.4  volts  in  either  the  positive  or  negative  direction 
The  calibration  levels  input  to  the  system  are  based  upon 
this  maximum;  signals  greater  than  1.4  volts  will  saturate 
the  amplifier  and  the  playback  will  simply  be  a  straight 
line  at  the  maxima.   Yet  on  the  recordings  for  virtually 
all  of  the  strain  gages  which  survived  to  produce  intel- 
ligible records,  the  output  was  greater  than  the  1.4  volt 
limit,  but  no  concurrent  indication  of  amplifier  saturation. 
Figures  6.2  and  6.3  are  offered  as  typical  representations 
of  this  phenomena  for  each  of  the  two  plates. 

A  number  of  possible  sources  for  the  strain  gage  error 
were  considered: 

1.  The  calibration  signals  were  improperly  set. 

2.  The  signals  were  being  altered  either  as  they  were 
being  input  to  or  output  from  the  Fourier  Analyzer. 

3.  The  playback  amplifiers  on  the  tape  recorders  were 
out  of  calibration. 
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Figure  6.2   Untranslated  Strain  Output  for  Tee-Stiffened 
Plate  (Gage  11) . 
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Figure  6.3   Untranslated  Strain  Output,  Rectangular- 
Stiffened  Plate  (Gage  15). 
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4.   The  calibration  signal  was  being  -input  to  the  signal 
conditioners  at  the  wrong  location. 


'  o 


The  first  three  of  these  postulated  sources  of  error  were 
readily  discarded.   The  calibration  signals  recorded  with 
the  strain  gage  signals  were  precisely  those  set  into  the 
system  prior  to  each  shot;  thus  ruling  out  the  possibility 
that  either  the  calibration  levels  were  incorrect  or  that 
the  playback  amplifiers  were  out  of  calibration.   To  check 
the  second  point,  a  series  of  calibrated  sine  waves  were 
input  to  the  Fourier  Analyzer  using  the  same  connections  and 
equipment  settings  used  for  the  strain  gage  data  recordings. 
These  calibrated  signals  were  recorded  on  disk  and  then 
output  to  the  plotter  in  exactly  the  same  manner  as  the 
strain  gage  signals  had  been.   In  no  case  did  the  Fourier 
Analyzer  alter  the  calibration  of  the  signal. 

Whether  the  remaining  suggested  cause  was,  indeed,  the 
source  of  error  or  whether  the  actual  cause  is  something 
entirely  different  cannot  be  determined  at  this  point. 
Whatever  the  actual  cause,  the  effect  remains  the  same;  the 
strain  magnitudes  derived  from  the  two  underwater  tests 
cannot  be  used  for  comparison  with  the  results  predicted  by 
the  USA-STAGS  code. 

Although  the  magnitudes  of  the  test  derived  strain  data 
are  clearly  unuseable ,  the  general  shape  of  the  strain 
curves  should  be  essentially  unaffected.   In  an  attempt  to 
resurrect  the  comparative  nature  of  this  study,  it  was 
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determined  that  the  shape  of  the  curves  from  the  tests 
should  be  contrasted  with  the  strain  curves  predicted  by  the 
USA-STAGS  code.   This  attempt,  too,  was  doomed  to 
f rust  rat  ion . 

Figures  6.4,  6.6,  and  6.8  are  representative  samples  of 
the  strain  data  calculated  by  USA-STAGS  functioning  as  a 

post-shot  analyzer  of  the  tee-stiffened  plate.   Figures  6.5, 

49 
6.7,  and  6.9  are  the  corresponding  strain  gage  outputs. 

Similarly,  Figures  6.10,  6.12,  and  6.14  are  USA-STAGS 
outputs  for  the  rectangular-stiffened  plate;  and  Figures 
6.11,  6.13,  and  6.15  are  their  corresponding  strain  gage 
histories . 

While  similarities  exist  at  some  points  between  the 

actual  strain  histories  and  the  USA-STAGS  predicted  results, 

50 
in  general  there  is  very  little  correlation.     This  -was 

unexpected  since  earlier  discussions  with  personnel  at 
Lockheed  had  indicated  that  extremely  high  degrees  of  corre- 
lation had  been  obtained  when  the  code  results  had  been 
compared  to  with  tests  conducted  by  the  Defense  Nuclear 
Agency  and  others. 


49 

Unfortunately,  after  all  of  the  problems  encountered 

with  the  strain  gages,  there  were  not  nearly  as  many  of 

these  as  might  have  been  desired. 

50 

Once  again  it  should  be  emphasized  that  the  shapes 

alone  of  the  curves  can  be  compared.   As  already  explained, 

the  magnitudes  of  the  strain  gage  records  are  not  correct. 
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Figure  6.4   Tee-Stiffened  Plate,  Gage  2  (USA-STAGS) 
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Figure  6.6      Tee-Stiffened  Plate,  Gage  8  (USA-STAGS) 
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Figure  6.7   Tee-Stiffened  Plate,  Gage  8 
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Figure  6.8   Tee-Stiffened  Plate,  Gage  11  (USA-STAGS) 
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Figure  6.9   Tee-Stiffened  Plate,  Gage  11 
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Figure  6.10   Rectangular-Stiffened  Plate  Gage  1  (USA-STAGS) 
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Figure  6.11   Rectangular-Stiffened  Plate,  Gage  1 
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Figure  6.12   Rectangular-Stiffened  Plate,  Gage  15  (USA-STAGS) 
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Figure  6.14   Rectangular-Stiffened  Plate,  Gage  18  (USA-STAGS) 
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Any  number  of  possible  solutions  exist  for  this  perform- 
ance, of  course;  including  the  extremes  of  either  USA-STAGS 
or  the  strain  histories  being  totally  wrong.   In  this  case, 
however,  utilizing  the  extended  capabilities  of  the  code 
itself  may  have  provided  the  most  plausible  solution. 

The  USA  postprocessor,  POSTPR ,  has  the  capability  of 
calculating  the  fluid  pressures  at  the  fluid  structural 

interface.   These  results  may  then  be  displayed  in  either 

51 
tabular  or  graphical  form/    When  this  was  done  for  the 

cases  studied  here,  a  very  interesting  result  was  obtained. 

Figures  6.16  and  6.17  represent  the  predicted  fluid 
pressures  at  the  center  of  the  tee-stiffened  and  rectangular- 
stiffened  plated,  respectively.   As  can  be  seen,  USA-STAGS 
predicts  that  intense  negat  ive  pressures  will  be  generated 
within  a  very  short  time  after  the  arrival  of  the  shock 

wave.   That  is,  USA-STAGS  is  predicting  massive  local 

52 

cavitation  across  the  faces  of  these  plates. 

Viewed  in  that  light,  the  reason  for  the  seemingly  poor 
performance  of  USA-STAGS  in  predicting  the  strain  responses 


51 

As  noted  earlier,  no  graphical  capability  exists  on 

the  VAX/VMS  system  at  the  Naval  Postgraduate  School.   These 

post-shot  analysis  problems  were  run  at  Lockheed  Research, 

Palo  Alto.   It  was  they  who  most  generously  provided  the 

accompanying  plots. 

52 

In  actuality,  local  cavitation  would  be  marked  by  near 

zero  pressures  as  the  water  flashed  into  steam  to  fill  the 
void  created  by  the  sudden  displacement  of  the  plate.   Since 
USA-STAGS  assumes  a  homogeneous  medium  this  must  necessarily 
be  represented  by  negative  pressures. 
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Figure  6.16   POSTPR  Predicated  Fluid  Pressures  for  Tee-Plate 
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can  be  better  understood,  and  even  the  strain  gage  responses 
themselves  become  more  understandable. 

For  example,  in  Figure  6.11  an  initial  compressive  strain 
can  be  observed  as  the  plate  reacts  to  the  incident  pres- 
sure.  This  does  not  occur  immediately  but  is  delayed  some- 
what as  the  inertia  of  the  plate  is  overcome.   In  this 
instance,  the  initial  compressive  strain  is  mimicked  almost 
exactly  by  the  response  predicted  by  USA-STAGS.   The  strains 
then  return  to  zero  or  become  tensile  as  cavitation  occurs. 
As  would  be  expected,  the  responses  predicted  by  the  code 
are  much  more  violent  then  those  seen  in  the  strain  gage 
histories.   This  is  so  because  they  are  activated  by  large 
values  of  assumed  pressure  rather  than  the  actual  zero  pres- 
sure that  would  occur  with  cavitation. 

The  other  plots  seem  to  confirm  this  basic  trend: 
Initial  zero  or  compressive  strain  depending  upon  the  loca- 
tion of  the  element /strain  gage  followed  by  a  positive  trend 

53 

in  the  strain  history  as  cavitation  develops. 

The  actual  mechanisms  involved  in  the  results  that  seem 
to  be  indicated  here  are  not  fully  understood.   The  implica- 
tion, however,  is  that  the  velocity  of  the  plate  as  it  is 
being  deformed  is  much  greater  than  that  of  the  fluid.   A 
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It  is  suggested  that  Figures  6.9  and  6.15  also  adhere 

to  this  basic  pattern,  but  are  subject  to  the  polarity 

problem  mentioned  earlier.   It  is  believed  that  these  plots 

should  actually  be  inverted. 
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discontinuity  is  therefore  being  created  in  the  fluid/ 
structural  interface.   This  discontinuity  is  manifested  by 
the  pressure  drop  which  causes  the  local  cavitation. 

The  strain  gage  histories  seem  to  imply  that  local  cavi- 
tation is  indeed  occurring.   Whether  that  is  in  fact  the 
case,  or  whether  the  cavitation  predicted  by  USA-STAGS  is 
simply  a  function  of  the  way  that  the  plates  are  modelled  is 
something  that  must  still  be  determined,  however. 

In  these  studies,  the  plate  alone  has  been  modelled. 
The  structure  that  is  "known"  to  the  computer  code  can 
therefore  be  conceived  as  infinite  in  extent  and  must  neces- 
sarily bear  the  full  brunt  of  the  pressure  wave.   In  actu- 
ality, this  is  not  so.   The  test  plate  is  finite  in  size  and 
pressure  relief  will  occur  around  its  boundaries.   Since  the 
breadth  of  the  shock  front  is  on  the  order  of  five  feet  for 
charges  of  the  weight  used  in  these  tests  (this  is  a  func- 
tion of  the  speed  of  sound  in  the  fluid  and  the  pressure 
decay,  it  will  vary  with  charge  weight),  the  entire  assembly 
will  be  enveloped  by  the  pressure  front.   This  means  that, 
almost  instantaneously  after  the  incident  pressure  is  felt 
by  the  test  plate,  the  pressure  will  begin  to  relieve  around 
the  edge  of  the  plate  and  a  negat  ive  pressure  will  be 

experienced  on  the  backs  of  the  flanges,  and,  soon  after,  on 

54 
the  rear  of  the  backing  structure.     Therefore  the  entire 
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Complex  interference  patterns  can  also  be  expected 

137 


force  of  the  incident  pressure  wave  may  not ,  in  fact  ,  be 
impinging  upon  the  test  plate;  since,  while  that  wave  is 
"pushing"  on  the  front  side  of  the  structure,  it  is  also 
"pulling"  on  the  back. 

Whether  this  decrease  in  force  would  be  sufficient  to 
eliminate  cavitation  is  an  open  question.   More  testing  must 
be  done.   For  the  USA-STAGS  code  to  be  able  to  predict  this, 
however,  it  is  not  enough  to  model  only  the  test  panel.   The 
entire  assembly  of  test  panel,  milled  edge,  and  backing 
structure  must  be  included.   If  local  cavitation  is  still 
predicted,  then  additional  action  will  be  required  because 
Doubly  Asymptotic  Approximation  which  is  the  basis  for 
the  USA  code  is  not  valid  where  discontinuities  such  as 

cavitation  exist  between  the  fluid  and  the  structural 

55 

surface    [Ref .  17] . 

C.   CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  STUDY 

Clearly,  problems  encountered  both  in  conducting  under- 
water testing  and  in  using  the  USA-STAGS  code  have  precluded 
any  firm  conclusions  from  being  made  about  the  predictive 
capabilities  of  the  code.   Thus,  the  primary  objective  of 
this  study  has  not  been  met.   It  is  in  the  fulfillment  of 
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Lockheed  has  developed  a  processor  ,  USA-STAGS  CFA 

(Cavitated  Fluid  Analyzer)  which  can  be  used  to  model  the 

cavitated  region.   The  process  that  must  be  followed  is  to 

model  layers  of  fluid  at  increasing  intervals  from  the  plate 

until  no  further  cavitation  has  been  found.   At  that  point, 

the  standard  DAA  approximation  may  be  applied. 
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the  second  objective  presented  in  the  introduction,  that  of 
providing  guidelines  for  future  study,  thai;  definite  contri- 
butions can  be  made. 

In  the  realm  of  underwater  explosion  testing,  some 
specific  problem  areas  were  noted  in  the  preceding  chapters. 
Solutions  to  some  of  these  (e.g.,  strain  gage  choice  and 
attachment)  were  presented  in  the  text  of  those  chapters. 
Others  of  these  problems  are  of  a  more  fundamental  nature 
and  must  be  discussed  separately. 

First,  a  satisfactory  plan  of  testing  cannot  be  devel- 
oped without  an  adequate  source  of  appropriately  sized 
explosive  charges.   The  inability  to  match  charge  weight  to 
test  platform  dimensions  makes  it  impossible  to  use  any 
computer  code  in  it-s  preshot  capacity  and  impacts  negatively 
upon  the  reproduceabil it y  of  test  results.   Either  the 
charge  weight  must  be  standardized  or  the  test  platform 
dimensions  must  be  altered  to  accommodate  changes  in  charge 
weight.   In  either  case,  the  most  apparent  solution  seems  to 
be  a  carefully  orchestrated  test  plan  with  long  lead  times 
before  each  test  shot.   If  a  match  is  not  possible,  it  is 
suggested  that  the  shot  should  probably  be  aborted.   That 
valuable  information  was  obtained  from  such  mismatched 
conditions  in  this  study  is  considered  to  be  more  a  charac- 
teristic of  a  program  in  its  developmental  stages  than  of 
well-conceived  experimental  practice.   Until  more  experience 
is  gained,  it  appears  that  the  rationale  dictating  use  of 
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the  eight  pound  TXT  charge  as  the  standard  [Ref.  3]  is  still 
valid.   It  is  therefore  recommended  that  that  weight  charge 
be  maintained  if  at  all  possible. 

Second,  the  lack  of  adequate  strain  gage  data  has  been 
an  incapacitating  flaw  in  the  test  results.   Whether  the 
poor  quality  of  that  data  has  been  the  result  of  errors  in 
the  equipment  set-up  at  the  West  Coast  Shock  Test  Facility, 
the  interface  between  that  equipment  and  the  strain  gages  on 
the  test  plate,  or  some  other  cause  is  uncertain.   What  i_s_ 
certain  is  that  the  present  arrangement  of  sending  the  test 
plate  and  backing  structure  to  WCSF  a  few  days  before  the 
test  shot  for  attachment  and  then  allowing  only  a  few  hours 
immediately  before  the  shot  for  equipment  testing  is 
insufficient.   Development  of  an  in-house  instrumentation 
capability  at  the  Naval  Postgraduate  School  is  considered  to 
be  essential.   This  would  permit  the  gages  and  electronic 
equipment  to  be  integrated  into  a  single  system.   Thorough 
testing  could  be  conducted  before  the  actual  date  of  the 
test;  and,  equally  important,  responsibility  for  the  outcome 
of  the  test  would  reside  in  a  single  organization. 

Changes  in  the  electronic  instrumentation  have  already 
been  recommended  [Ref.  3].   Enactment  of  those  changes  would 
undoubtedly  improve  the  data  responses.   A  few  other  modifi- 
cations suggested  by  the  results  of  this  study  are  consid- 
ered equally  important,  however. 


140 


1.  Rather  than  using  calibration  resistors,  the  strain 
gage  bridge  should  be  balanced  using  strain  gages 
identical  to  those  on  the  test  plate.   These  should 
be  attached  to  an  aluminum  plate  in  the  terminal  box. 
This  has  two  effects:   First,  the  balancing  gages 
would  be  subjected  to  the  same  changes  in  temperature 
as  the  plate  gages.   Second,  the  balancing  leg  of 

the  bridge  would  be  of  the  same  length  as  the  measure- 
ment leg.   Combined,  these  would  make  the  bridge 
far  less  sensitive  to  environmental  changes.  '  A 
zero  balance  could  be  maintained  and  calibration 
fluctuations  would  be  less  likely. 

2.  A  zero-time  triggering  signal  should  be  superimposed 
upon  all  instrument  channels.   This  would  permit  an 
exact  assessment  of  the  time  of  initial  plate  response 
to  the  incident  wave. 

3.  Further  shielding  and  grounding  of  the  pressure  gage 
cables  should  be  attempted  to  prevent  the  saturation  of 
the  pressure  gage  amplifiers  through  currents  induced 
by  charge  detonation. 

From  the  results  that  have  been  explored  in  the 
preceding  section  of  this  chapter  it  is  clear  that  certain 
questions  need  to  be  raised  about  the  validity  of  the  flat 
plate  model  that  has  been  used  here.   These  same  questions 
should  be  asked  about  the  ability  of  USA-STAGS  or  any  other 
code  to  handle  that  model  under  the  conditions  imposed. 
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In  particular  the  question  of  whether  the  local  cavita- 
tion predicted  by  USA-STAGS  is  actually  occurring  or  is 
simply  a  result  of  the  choice  of  model  must  be  addressed. 

It  is  recommended  that  this  study  should  be  redone,  but 
with  the  backing  structure  included  as  a  part  of  the  basic 
model.   POSTPR  could  then  be  run  to  determine  whether  cavi- 
tation was  still  occurring  across  the  plate.   If  cavitation 
was  still  being  indicated,  then  either  a  new  model  would 
have  to  be  found  or  the  cavitated  region  would  have  to  be 
modelled  using  a  code  capable  of  handling  it. 

Whatever  the  findings  of  continued  tests,  it  may  have  to 
be  recognized  that  increasing  the  scope  of  the  models 
studied  may  put  USA-STAGS  beyond  the  reach  of  masters-level 
research  at  the  Naval  Postgraduate  School.   USA-STAGS  is  a 
diverse,  highly  capable  code  which  has  the  reputation  for 
high  levels  of  accuracy.   The  price  that  is  paid  for  such 
capability,  however,  is  size  and  complexity.   Given  the 
limited  time  and  computer  resources  available  to  most 
students  at  NPS ,  USA-STAGS  may  be  beyond  their  capacity.   If 
such  research  is  undertaken,  it  is  most  emphatically  recom- 
mended that  a  large  block  of  time,  weeks,  perhaps,  be  spent 
at  Lockheed  Research  in  Palo  Alto,  California  learning  the 
code  and  perfecting  models. 


142 


APPENDIX  A 
SOURCE  DATA 

This  appendix  contains  the  information  listings  supplied 
with  the  STAGS  and  USA  source  tapes.   Section  A  is  a  full 
listing  of  the  contents  of  the  STAGS84  tape  along  with 
instructions  for  mounting  the  tape,  command  files  for 
creating  and  linking  executables,  and  sample  command  files 
for  running  the  code.   Section  B  contains  similar  listings 
for  the  USA  source  tape  excepting  that  no  command  files  are 
listed.   USA  places  the  link  commands  on  the  tape  itself  for 
immediate  access.   Command  files  for  running  the  modules  are 
identical  in  construction  to  those  for  STAGS  and  can  be 
adapted  from  those  in  Section  A. 
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A.   THE  STAGS84  SOURCE  TAPE 


»Zi-.  SIGN   Au^uST  l9ti« 
STASSC1  PROGRAM  T~R£ 
VAX /VMS  OPERATING  SYSTEM 


TAPE  FORMAT 


9-Track 
1600  BPI 
Labeled  ( STAGS4 ) 


TAPE  CONTENTS 

File  Name  Type 

1  INFO  VAX  Tex  t 

2  STAGSU  FOR  Source 

3  3LICMF  MAR  Source 

4  J0E3TA. MAR  Source 

5  STAGS1  FOR  Source 

6  MOUNTEL. FOR  Source 
6  STAG32  FOR  Source 
9  STAPL.  FOR  Source 
9  POSTP  FOR  Source 

10  STAGS1  0L3  Liorary 

11  STAGS2  0L3  Library 

12  STAGSU  0L3  Library 

13  POSTP.  OLB  Library 

14  STAPL.  OLB  Library 

15  NICE  OLB  Library 

16  TP  FOR  Source 

17  TP.  OBJ  Ob  ject 

18  ZSYS.  FOR  S'ource 

19  RES.  FOR  Source 

20  RES  OBJ  Object 

21  ZSYS. OLB  Library 
22-36  CASE*.  *.  •  List 


Description 

Table  of  contar.  rsi   procedure-; 
FORTRAN  prog-ams  for  STAGS  utility  library 
MACRO  program  far  STAGS  utility   1  i  b  r  a  r  y 
MACRO  program  :or  STAGS  utility   liDrary 
FORTRAN  prog-am?  forSTAGSl   library 
FORTRAN  programs  'or  the  MOUNT  elements 
FORTRAN  programs  for  5TAGS2  library 
FORTRAN  prog-am?  for  STAPL  library 
FORTRAN  prog-am?  for  POSTF  library. 
Library  of  S  ~  —  C-  5  1  programs 
Library  of  5T-GS2  programs 
Library  of  S_i-GE  utility  routines. 
LiDrary  of  POSTP  programs 
Library  or  STAPL  plot  programs. 

Thurston  Processor. 

Thurston  Processor. 

Z-system  utilities  (for  TP ) . 

Reduced  Equation  Solver 

Reduced  Equation  Solver 

Z-system  utilities  object  library 

Input  and  output  for  5  sample  cases 


3.     LOADINC  FROM  TAPE 

Us«  the  following  control  cards  to  copy  all  files  to  disc. 

*  MOUNT  MTAO:  STAGS4  TAPE 

*  COPY  TAPE:  *.  *.  *  C] 

Files  may  also  be  read  indi  dually  from  tape  ujitn  the 
following  control  cards  (e  g.  ,  to  load  STAGS!  0L3  '>  . 

%    COPY  TAPE- STAGS1. OLB  C] 

Complete  sample  case  output  may  be  oo  "aired  as  follows: 

*  PRINT/HEAD  CASE«  *  * 
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4  PROCEDURES 

A    LiBR  COM    (Compile  programs  an;  c  r  s  ■=  "  s  1  1  t  ^  :  i  ;  s  i 

i    "OF  3TAG5U  FOR /NOLI ST 

*  LIBRARY/CREATE  STAGSu  QL3  STAGSU 
$  MACRO  3LIGMF  MAR 

*  MACRO  JOSSTA  MAR 

*  LIB  STAGSU  SLIOMF.  JQ3E-TA 

5  POP  STAGS1  .FOR/NOLIST 

5  LIBRARY/CREATE  STAGS  1  GL3  5TAGS1 

s  FOR  MQUNTEL  FOR/NOLI 3T 

*  LIB  STAGS1  MOUNTEL 

t  FOP  3TAGS2  FOR/NQLIST 

S  LIBRARY/CREATE  STAGS2  OLB  STAG32 

S  FOR  MOUNTEL.  FOR/NQLIST 

S  LIB  STAGS2  MOUNTEL 

*  FOR  5TAPL  FOR/NOLIST 

*  LIBRARY/CREATE  STAPL  OLB  STaPL 

4  ^GR  POSTP  FOR/NOLIST 

*  LIBRARY/CREATE  POSTP  OLB  POSTP 

»  FOR  TP  FOR/NOLIST 
t    FOR  ZSYS  F0P/N0LI3T 

*  LIB/CREATE  ZSYS  ZSYS 

*  FOR  RES. FOR/NOLIST 

B.  LSTAGS1.COM    (Load  STAGS1) 

*  L INK /MAP=STAGS1 /BRIEF /EXEC UTAELE=STAGS1- 

STAGS1  0LB/LI3/INCLUDE=MAIN1, - 
STAGSU. OLB/LIB 

C.  LSTAGS2.COM    (Load  STAGS2) 

5  LINK/MAP=STAGS2/ BRIEF /EXECUT A3LE=STAGS2- 

STAG32.  0L3/LIB/INCLUDE=MAIN2,  - 
STAGSU  OLB/LIB 

D.  LSTAPL.COM    (Load  STAPL) 

*  LINK/MAP=STAPL/BRIEF/EXECUTABLE=STAPL- 

STAPL.  OLB/LIB/INCLUDE=STAPL, - 

STAGSU.  OLB/LIB. - 

(Fil«  containing  Local  plot  routines) 

E.  LPOSTP  COM    (Load  POSTP ) 

$  LI NK/MAP=POSTP/ GRIEF /EXECUT ABLE=POSTP- 

POSTP  OLB/LIB/  1NCL,JDE=P0STP,  - 
5TAGS1  OLB/LIB, STAGS2  CLB/LIB-- 
STAGSU  GL5/LIB, NICE.  QL3/LIB 


145 


STAG'S!  COM 


>=:: i     - . ijii 


s  ^aSISN  CASE  I  INP  F OR 00 5 

5  ASSIGN  Ca£E1  OUT  FQROOc 

s  ASSIGN  CA5E1  MOD  F0R021 

s  RUN  ETAGS1  EXE 

4  PRINT  CASE!  OUT 


g.   stags  con 


(Execute  STAGS1  anc  STAGS2) 


*  ASSIGN  CA5E1  INF  F0P005 
$  ASSIGN  CASE1.  OUT  FORCCS 

*  ASSIGN  CASE!  MOD  FOR 021 
$  PUN  STAGS  1  EXE 

*  ASSIGN  CASE1  RST  F0RG2G 

*  ASSIGN  CA5E1.S0D  F0RC22 

*  PUN  STAG32  EXE 

*  PRINT  CASE!  OUT 


(optional  if    model  plots  o  e  s  l  r  e  d  ) 

(previous  solution  data,   i-e^d   tor  restart) 
(optional  if  soiuuon  data  s  n  o  u  i  d  oe  raved) 


H 


POSTP  COM    (Execute  STAGS1  ana  PQSTp ; 


*  ASSIGN  CASE4  INF  F0RGG5 

*  ASSIGN  CA5E4  OUT  F0POG6 
s  RUN  STAGS  1.  EXE 

*  ASSIGN  CA3E4  PPN  F0R005 
i  ASSIGN  CASE*  POP  F0R006 
J  ASSIGN  CA5E4  RST  F0R020 
$  ASSIGN  CASE*  SOD  F0R022 

*  RUN  POSTP  EXE 

*  DEL  FOR009  OAT 


I.   STAPL  COM    (Execute  STAPD 

*  ASSIGN  CASE1  PIN  F0R005 
S  ASSIGN  CASE1  RST  FOF020 
5  ASSIGN  CASEl  MOD  F0R021 

*  ASSIGN  CASEl. OUT  FQR006 

*  RUN  STAPL. EXE 

*  (Execute  plot  file  on  local  plot  device 


(May  be  omitted  if    FORGO?.  DAT  saved; 
(previous  solution  data  from  CASE*  SOD) 
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B.   THE  USA  SOURCE  TAPE 


'-.•'MDERwATslP  5HQCK  '-NA_':-.  .'5  (  US<- 
VAa  TaPE  LOADING  INSTRUCTIONS 

'f  1:      tape      C  5  n  t  a  1  P  s      Symbolical       ir.r  iuCJ      .;  •   0  C  e  0  'J  T  ■?  S  •       "C  j;ct      I  1  -j  r-  g  r  l  > 

»  '  5  c  u  1 5  c  1  ?  5  .  ana  manuals  -or  a  number  at  US  .-•  :ro:;riors  All  f  i  i  e  ;  r-  a  v  ■ 
Oe  =  n  :r»at?(  with  COPY  c  orr.ma  n  a  3  ana  tan  De  "strisvsd  in  the  same  un  n<?; 
with  only  a  'e«i  lines  of  input  a*  the  t  e  r  m  i "  -j  1  Examination  of  the  tap' 
directory  will  reveal  t  n  a  t  tn?  files  njve  .  a  -  i  a  u  s  c  y  c  i  e  n  u  m  o  e  r  s  inj  t  n  : 
13  tne  [»y  to  Keeping  tnose  files  together  - ,-,  at  Belong  to  a  particular 
processor.   The  following  taoie  contains  t  K  l  s  identification  coae 

USA  PROCESSOR  CYCLE 


FLUMAS  1 

AUGMAT  2 

TIM  I  NT  2 

post=>p  4 

MODLIB  5 

To  load  the  TIM  INT  files,   for  exam,  pie.  50  ~:  :n»  directory  'jn?r=  -.hey 
are  to  reside  and  simoly  do  a  COPY  TAPE  *  ■'*.  3  C]  command  .^riere  TAPE  is 
the  internal  name  that  has  aeen  assigned  t;  the  magnetic  tape  unit   This 
process  is  simply  repeated  until  all  a  e  s  1  r  e  c  cirectories  r,  a  » ■■  e  Peen 
loaded 

The  3  character   extension'  on  the  file  names  denotes  tne  tope  of 
file  according  to  the  following  table 

EXTENSION  TYPE 


. PDP  INCLUDE  PROCEDURES 

. INC  INCLUDE  PROCEDURES 

.  VAX  SYMBOLICS 

. OLB  OBJECT  LIBRARIES 

COM  LINK  STATEMENTS 

.  EXE  EXECUTAQLES 

.  UAL  MANUALS 

For  example,  AUGMAN  UAL  is  a  file  containing  the  AUGMAT  manual  while 
LINKUSAS.COM  is  the  LINK  statement  for  USA-STAGS.   It  should  De  noted 
that  all  LINK  statements  must  be  modified  slightly  as  the  directories 
referenced  in  them  reflect  names  that  are  in  use  here  at  Lockheed  and 
contain  users  names  that  will  not  exist  at  other  VAX  installations 
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U  1  r  s  r  ■ 


ADn-SS 

EXE- 

1 

AMD«TA 

POP, 

1 

ArtGEOl. 

FDP; 

1 

AMGE02 

FDP, 

1 

AMGE03 

PDP, 

1 

AMGE04 

PDP, 

1 

FLULIB 

OLB- 

1 

FLUMAN 

UAL; 

1 

FLUMAS 

FOR, 

1 

LINKFLU  CQf 

,  1 

IMAMFLU 

PDP. 

1 

AMGEC1 

PDP. 

-i 

c 

AMGE03 

FDP, 

-1 

AIJGLIB 

OLD, 

2 

AUOMAN. 

UAL, 

AUGMAT 

FOR, 

FIi_DAT 

PDP, 

e. 

LINK.AUG 

!  COM,  2 

NAMAUG 

PDP. 

- 

c 

PPE^RQC 

EXE.  2 

TPAI43F 

PDP. 

3 

& 

C KG DAT 

PDP; 

3 

COMADD. 

PDP; 

-3 

CSHPAC. 

INC; 

3 

DATFIL 

PDP, 

3 

DATINT. 

PDP; 

3 

INTLIB 

OLB. 

3 

INTMAN 

UAL. 

3 

LIN* USA  COM,  3 

NAMINT 

PDP; 

3 

TIMINT 

FOR, 

3 

UNUASH 

EXE, 

3 

COMADD. 

PDP, 

4 

DATFIL. 

PDP 

4 

LINKPOS.  COM;  4 

NAMINT 

PDP 

4 

PQSLIB 

OLB 

4 

PQ3MAN 

UAL 

4 

P03PP0 

FOR 

4 

POSTPR 

EXE 

4 

CHGDAT 

PDP 

5 

COMADD 

PDP 

5 

DATFIL 

PDP 

5 

DATINT 

PDP 

5 

LINKUSAS. COM; 

MQDLIB 

FOR 

5 

MQDLIB 

OLB 

5 

NAMINT 

PDP 

5 

USA'S  EXE-  5 

US TL  I  2 

OLB 

5 

UST3YM 

FOR 

5 

SI03UP 

MAR 

8 

CC  FOR 

g 

CGMERR 

PRC 

3 

DMGASP 

FOR 

8 

DMGASP 

ODJ 

5 

n  /  1    _i  ,- 

(:■-:£-- 

1  wr  • 

uu 

1/  1 

17-hA'- 

1  934 

'„':- 

00 

1/  1 

4- JUN  — 

1984 

r  j* 

00 

1/  1 

4-JUN- 

193- 

1 ; 

00 

1/1 

19-JUL- 

1  98  '. 

u  ■ 

00 

1/1 

4— JUN- 

1934 

.   ■„ 

00 

i9/69 

6-3EP- 

1984 

00 

00 

3  7  /  3  7 

6-5EF- 

1  jg^ 

00 

00 

74/74 

6-3EP- 

1  93- 

''jC  - 

00 

1/1 

6-3EP- 

193, 

M*  ' 

00 

1/1 

22-APP- 

1931 

>■-■ 

00 

1/1 

4-JUN- 

1934 

oc 

00 

1/1 

22-JUL- 

198". 

oc . 

00 

25/25 

6-SEP- 

1984 

vj -,-  '. 

00 

12/12 

13-MAY- 

i9s: 

00  ' 

00 

23/23 

fc-SEP- 

1934 

-■G 

00 

1/1 

23-MAY- 

1934 

00 

00 

1/1 

6-3EP- 

1984 

00 

00 

1/1 

22-APR- 

1981 

r,Q 

00 

42/42 

6-SEF- 

198  4 

r.  f, 

00 

1/1 

29-MAY- 

1984 

'-"  -' 

00 

1/  1 

30-MAR- 

1984 

00 

00 

1/1 

30-MAF- 

1984 

00' 

00 

1/1 

29-MAP- 

198" 

00 

00 

1/  1 

15-JUN- 

1  =5S- 

00 

00 

1/  1  . 

30-MAR- 

1984 

oc 

00 

53/53 

6-SEP- 

198-1 

0  0 

00 

22/22 

3-N0V- 

1983 

00 

00 

1/1 

6-SEP- 

1984 

00 

00 

1/1 

30-MAR- 

1984 

00 

00 

58/58 

6-SEP- 

1984 

00 

00 

70/70 

6-SEP- 

1984 

00 

00 

1/1 

30-MAR- 

1984 

00 

00 

1/1 

30-MAR- 

1984 

00 

00 

1/1 

30-MAR- 

1984 

00 

00 

1/1 

30-MAR- 

1984 

00 

00 

38/38 

30-MAR- 

•1984 

00 

00 

15/15 

27-MAR- 

1934 

00 

00 

41/41 

29-MAR- 

•1934 

00 

00 

61/61 

30-MAfi- 

1984 

00 

00 

1/1 

5-APR- 

1983 

00 

00 

1/1 

18-MAR- 

•1983 

00 

00 

1/1 

15-JUN- 

•1934 

00 

00 

1/1 

15-JUL- 

19e2 

oc 

00 

1/1 

6-SEP- 

•1984 

00 

00 

10/10 

6-5EP- 

1984 

00 

00 

9/9 

6-SEP- 

1984 

00 

00 

1/1 

22-APR- 

•1982 

00 

00 

122/223 

6-SEP- 

■la8-i 

00 

00 

7/7 

6-SEP- 

1984 

CtO 

00 

1  1  /  1  1 

6-5EP- 

-1984 

00 

00 

9/9 

15-AUG- 

•1981 

00 

00 

1/1 

9-5EP- 

■1982 

00 

00 

1/1 

6-JAN- 

■1932 

00 

00 

18/18 

8-JAfJ- 

la82 

00 

00 

12/12 

6  —  J  A  V*  - 

■1982 

00 

00 
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I.  J  f '!  L  jr       r   R  l"_  .  2 

1/  1 

t  -  .  - '  -  - 

i  -  - .. 

_■  . 

■J'J 

ig-'.om  -?:.  5 

1  '  1 

O— JAW- 

■  —  —  r 

_ 

Z  - 

[  G  M  D  A  T  P  R  C  •  5 

1,  1 

s- JAN- 

i  -  3  •: 

00 

I OH G c F  r  R  C .  3 

1/1 

~~. AN- 

'.951 

'   * 

J'Z 

rCHEQvi  =RC.  3 

i  •'  1 

G-jAi'J- 

■  ;c; 

.  . 

NICE  CLE ; 3 

lz.5/2-5 

10-Aw'G- 

195- 

.*.  ." 

00 

OUCH  TEE,  3 

1/  1 

21-NGV- 

1  95  j 

'■_■  . 

00 

3YM3IG  =RC; 3 

1/1 

6-JAN- 

1  931 

"  .". 

00 

5*l"iEuL  PRC;  3     "~- " 

6--JAN- 

1932 

lJ'-J 

00 

sYMDAT  PRC  ;  5 

i  n 

6-JAN- 

1932 

G'J 

00 

TE3TGA3P  TES. 3 

1/1 

6-JAN- 

1982 

00 

00 

Si*  "LINE  FOR;  9 

12/12 

2-AUG- 

19S2 

'j  r_! 

00 

SALINE  QL3.  9 

9/9 

10-AUG- 

1952 

u  'SJ 

Oo 

CYL  DA  A,  1  1 

3/3 

30-MAR- 

1984 

00 

CO 

CYL  FLU, 11 

3/3 

30-MAR- 

1934 

GO 

00 

CYL.  GEO,  1  1 

2/2 

30-MAR- 

1984 

.  j  ■"■ 

00 

CYL  P03. 1 1 

161/161 

l-APR- 

198- 

'j  0 

00 

CYL.  PRE,  1  1 

5/5 

30-MAR- 

L9S4 

0  0 

00 

CYL  RST, 1 1 

20/20 

1-APR- 

1984 

0'_ 

00 

CYL  USC, 1 1 

1/  1 

13-JAN- 

1934 

DC 

30 

CYL  USD. 1 1 

2/2 

13-JAN- 

193- 

■  ~>r. 

CO 

Ct'LAUG  INP;  11 

1/1 

17-FE3- 

1931 

00 

00 

C  VLAUG  OUT,  1 1 

30-MAR- 

198- 

GC 

00 

C  rL.PLU  INP,  11 

1/  1 

4-.JAN- 

193  4 

00 

CO 

CYLFLU  GUT, 1 1 

39/39 

30-MAR- 

193- 

00 

00 

CYL  I  NT.  INP,  1  1> 

1/  1 

1-APP- 

198- 

-,  J 

00 

CYLINT  out:  1 1  . 

80/30 

1-APR- 

1984 

00 

00 

CfLPQS  INF, 1 1 

1/1 

1-APP- 

198- 

GO 

00 

CYLPOS  OUT, 1 1 

113/113 

1-AFR- 

198- 

DC 

00 

CYLSTG  INP, 1 1 

1/1 

4-JAN- 

1934 

i"  r 

00 

CYL3TG  OUT,  1  1 

13/13 

1-APR- 

19S- 

00 

00 

P0R002  DAT, 1 1 

120/120 

1-APR- 

193- 

00 

CO 

FGRG22  DAT, 1 1 

15/15 

1-APP- 

1^3- 

(•>.;-". 

00 

BEAM  DAAi  15 

3/3 

3-JAN- 

1984 

00 

00 

EEam  GEO; 15 

1/1 

3-JAN- 

198-i 

CO 

00 

BEAM  PRE, 1 5 

2/2 

3-.JAN- 

1984 

00 

00 

BEAM. USC, 15 

1/1 

3-JAN- 

1984 

00 

00 

BEAM  USD; 15 

1/1 

3-JAN- 

1984 

00 

00 

BEAMAUG  INP;  15 

1/1 

3-JAN- 

1984 

00 

00 

BEAMFLU.  INP;  15 

1/1 

24-JAN- 

1984 

00 

00 

3EAMINE  INP, 15 

1/1 

3-JAN- 

1984 

00 

00 

BEAM  I MS  INP,  15 

1/1 

3-JAN- 

■1964 

00 

00 

3EAMPSE.  INP,  15 

1/1 

3-JAN- 

•1984 

00 

00 

EEAMPSS  INP,  15 

1/1 

3-JAN- 

■1984 

00 

GO 

BEAMSTG.  INP,  15 

1/1 

3-JAN- 

■1964 

00 

00 

Total  of  101  files-   1954/1954  blocks. 
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APPENDIX  B 
BASIC  SHELL  TYPES  FOR  USE  WITH  STAGS 

This  appendix  contains  a  detailed  listing  of  the  eleven 
basic  shell  types  used  with  the  STAGS  code.   Drawings  are 
included  to  show  numbering  conventions  for  edges  and 
corners.   This  information  can  be  particularly  useful  for 
discretizing  the  shell  unit  and  specifying  compatibilities 
wit h  .adjoining  shells.   (All  figures  from  Reference  13.) 
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RECTANGULAR  PLATE 


x.X 


ISHELL  =   3 


y.Y 


PROP(l)  =  X 

PROPU)    =    X 

PROP' 3)  =  Y 
PR0P(4)  =  Y^ 


Subscripts    refer  :o 
corner  r.umber. 


QUADRILATERAL  PLATE: 

t{   =  (1-XJ  •  (l-Y),    £,   =  X(1-Y).    f3  rXY,    f4  =X(1-Y] 

x(XY) 


y.Y 


E      £.  X. 


y(X  Y)  =    E     f.  Y. 
i.=  l    i     i 


PROP(l)  =  X 
PRCP(2)  =  Y 


PROPU)  =  X, 

paop(-i)  =  y; 


PROP(5)  =   X 


x,X 


PROP(o) 


3 


,1 


PROP(7)  =  X 
PROP(3)  -   Y~ 


ANNULAR  PLATE: 


2  2  2  2 

y~   t  ;~   -   r~   -  X 

x  =   t(X,    Y)   =   0 

y  =  g  (X,    Y)   =  X  sa  Y 

z  =  h  (X,    Y)  =  X  cos   Y 


PROP(l) 
FROP(2) 
PRCP(3) 
PR.OP(4) 


3   Rb 

-  a 

3  9b 


Liceral  displacement   posit: 
(  negative  x  ) '. 
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CYLINDER: 

2         2' 

y~  +  z     -  R~ 

x  =  f  (X.    Y)  =  X 
y  =  g  (X,    Y)  =  R 

z  =  h  (X,    Y)  =  R 


3  La 

COS 


V 

Y 


PROP(l) 
PROP(Z) 
PROP(3) 
PROP(4) 
PROP(5) 


=  X, 


Y, 


=   9. 


=  R 


L  =  X4  -  X1 


ISHZLI 


CONE: 

y2  *  z2  =  r2   =((X  -   X    ) 
tan  =   =  (Rh   -  Ral) 

x  =  £(X,    Y)   =  X  -   Xt 
y  =  g  (X,  Y)  «((X  -  Xt) 
z  =  h  (X,  Y)  =({X-  X    ) 


•    tan  a  -  R,) 
/  (x4  -  Xil 

-  R 


can  3 
tan  a 


a ' 

RJ 


CC3 


PROP(l)  = 

PROP(2)  = 

PROP(3)  = 

PROP(4)  = 

PROP(5)  = 

PROP(6)  - 


i! 

v 
R2 


R     may  be  >  Ru 


X,X       (Axis  of  revolution) 
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:?_E: 


2         2         2  2 

x    +  y     ■*■  z     =  R 

x  =  i(X,    Y)  =  -  R  =03  X 
y  =  g  (X,    Y)  =  R  sia  X  sia 
2  =  h  (X,    Y)  =  RsiaX  cos 


PROPfl)  =  X,    =  i 
1  a 

PROP'2)  =  X2  =  7b 
PROP(3)  =  Y,   =  g 
PRCP(4)  =  Y,  =  9^ 

4  0 

PROP(5)  =  R 


ISHTLJ-  =  3 


TORUS: 


x  =  £(X,    Y)  =   -Rq  cos   X 


y  =  g(X,    Y)  =  (R"   -  R,    3ia  X!   3ia  Y 
z  =  h(X,    Y)  =  (R,    fR.    sia  X)  c=a   Y 


PROP(l)  = 

Xl 

= 

7a 

PROP(2)  - 

X2 

= 

*b 

PROP(3)  = 

1  1 

= 

•5 
a. 

PRCP(4)   = 

Y, 

= 

5, 

PRO?(5)  =  R 
a 

PRCP(6)  r  R_ 
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[5HSLL,  =  3 


Z  L.LLP~T " 
■>  ■> 


co:;Z/C 


=   (x 


tan  3  =  (R  , 


*  R 


=     (X 


xL) 


za 

tan  y  ■  tan 


:a 

)/(X.-X, 


f(X,  Y)  =  X  -  Xi 
g(X.Y)   r  r     3  in  Ye 

h(X.  Y)  =  ry  cos   Ya 

z 

Y      =  ircua  (tan  Y  • 


*b 


PROP(l)  =  X. 

prop(2)  =  x: 

PROP(3)  =  Yi   =  9 
PROP(4)  =  Y, 
PROP(5)  =   R" 
PROP(6)  =  R£* 
PROP(7)  =  EUi 


R       mav  be  >  R    ,  . 
za  zb 


R      /R 

ya       za 


R      /R      ) 
za       ya 


ISKEL.L  =   10 


PARABOLOID: 


y     +•  z     =  r 
x  =  £(X,    Y) 


=  4R 


a 
Xi 


(x 


s> 


,yRa  (X  f  r  j 


z   =  2,yRa  (X  -  Rb)  •    cos 
PROP(l)  =  X, 


PROP( 
P.ROP(3 


X, 


I    "   -a 


PROP(4)  =   Y,   =  9b 
PROP(5)  =  Ra 
PROP(b)  =  Rs 


Axis  of  revolution) 
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1  I 


x  =  £(X,    Y)   =   -R   cos  Xc 
x 

y  =  g(X.  Y)  =  R   raiaX*  »La  Y 

z  =  h(X.   Y)   =  R       sia  Xa  cot  Y 


X     =  arctin  (taa  X  ■  ?.    /?.  ) 

y=     x 


PROP(4)  = 

Y-,  =  a. 
2       "b 

PROP(5)  = 

R 

X 

PROP(6)  = 

R 
7= 

HYPERBOLOID: 

7              ■> 

Y      '  Z 

(x  -  RJ" 

R- 

a 

R- 

0 

=  I 


x  =     X  -  X 


(R    /R  )  VRT 

a       o  d 


(X-X.-R    )■ 
1       c 


(R    /R  )  VR"  -  (X-X,  -R    )" 
a       o  d  1        c 


PROP(l)  =  X] 
PRQPC) = X, 
?R0P(3)  -  Y, 
PROP(4) 
PROP(5)  =  R 
PROP(6)  =  R 
PRCP(T)  =   R 


L    "    "a 
PROP(4)   =Y      i9 


155 


APPENDIX  C 
USA-STAGS  INPUT  FILES 

This  appendix  contains  the  input  files  used  for  the 
complete  USA-STAGS  runs  for  both  the  tee-stiffened  and 
rectangular-stiffened  plates.   These  are  listed  in  the  order 
in  which  they  would  normally  be  run: 

-  STAGS2  (Modal  Analysis)  Input  File 

-  STAGS1  Input  File 

-  Fluid  Mass  Input 

-  Augmented  Matrix  Input 

-  Time  Integrator  Input 

-  USA  postprocessor  (POSTPR)  Input 

The  input  files  for  the  tee-stiffened  plate  include  an 
example  of  STAPL  input.   The  files  for  the  rectangular- 
stiffened  plate  contain  an  example  of  POSTP  input  file. 
Since  these  last  files  are  readily  generalized,  it  was  not 
considered  necessary  to  provide  examples  for  each  case. 
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TEESTAG32.INP; 10 


9-3ED- 1 994  1 


3  2  4  4 
3  2  5  2 
1  0  S  I 
10. Eb  . 
11!$ 
.  1875 
1  1  3 
.125 
0  5  m 

6.  0 
0  0. 


STAGS2  OAT  A  FDR  TRANSVERSE  TEE  STIrFENER  (yOQAL  ANALYSIS) 
200000  100S3-1 

5  0  0  a  S  3-2 
10  2  5  3-3 
0.  5  C-l 
1  0  1  0000  $  0-2 
12  200.  0.  2  b  0  0.  3  0  - 3 
115  11  11  11  3  11  3  11  3  5  F  -  ! 
1  2  2  a  S  G-l-1  C0MPATA8ILI  TY  1-2 
1  2  3  a  S  G-l-2  C0MPATA3ILITY  1-3 
5  G-l-3  CO*-"  RATABILITY  3-4 
S  G-l-1  C0-MPATA8ILITY  3-5 
-1 

3  0.  .093  S  1-2  MATERIAL  PROPERTIES 
K-l-1  FLOOR  SHELLS 
5  K-2-1 

K-l-2  STIFFEN  £  R  SHELLS 
5  < -2-2 
-1-1 

.   3.  $  M-2A-1 
-.0938  J  M-5-1 
-1-1 

5  P-l-1  BOUNDS  1 
-1 
-1 

-1  -2 

.  9  .  $  M  -  2  A  -  2 
-.0933  S  M-5-2 
-1-2 

3  P-l-2  BOUNDS  2 
-2 
-2 

-1  -3 

.  .7313  3  M - 2  A  - 3 
3  .  .7813  5  M -  3  -  3 

0.  3  M-5-3 
-1-3 

5  P-l-3  BOUNDS  3 
-3 
-3 

-1  -a 

.  3.35?u  S  *-2K-a 
3  0.  3  .  3  5  9  u  .7313  I    m  -  3  -  4 

0  0  .  .0b25  5  M-5-U 
a  1  0  3  N  -  1  -  u 

3  2  5  5  P-i  -a  30'jnos  a 

S  3  -  1  -  4 
5  R-  1  -u 

2  3  M  -  1  -  5 
,  o.  2.6406  3.  5  M - 2  A  - 5 

1  0.  2  .  54U6  .7813  I  u-  3-5 
0  0.  .0525  3  M-5-5 

y  1 0  $  n  -  ■  1  -  5 

1  5  2  3  5  p-1-5  bounds  5 

0  3  Q  -  1  -  5 

0  5  R-l-5 


a  1 0 

4  6 
0  $ 

0  $ 
2  0 
0. 

1  0 
41  0 
u  2 


0 
0 
2 
0, 

1 

2  0 
410 

4  6 


3  N 

2  4 

a-i 

R-l 

S  H 

6.  3 
0. 

3  N 
2  6 

$  Q-l 
$  R-l 

2  S  M 
6.  0 

3  0. 
0  90. 

5  N 
2  6 
Q-l 
R-l 
$  M 
5.  3 
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Fiji! 
1   1   1  0  0  1  0 

0  0  4  4  3-2 


'R ANSWER SE  TEE  STIi 


:ZUi 


0. 

0 

0 

0  4  4  t- 

116  11 


1  *  B-l 


E-i 


STAGS  1  DA7 
6 
5 

10  2  5  3-3 
5  C-l 
1000  1  1  s  D-l 
2.  E-3  20  E-o  4 
:-2 

11  11  3  11  3  11  3  i  F-l 
5  G-l-1  COMPATAEILITY  1-2 

*  G-l-2  COMPATAEILITY  1-3 

*  G-l-3  COMPATAEILITY  3-4 
4  0-1-4  COMPATAEILITY  3-5 
1-1 

3  0.  .099  *  1-2  MATERIAL  PROPERTIES 


IPLAS" 


:o 


1  2  2 

1  2  3 

3  2  4 

3  2  5 

1  3  1 


10.  E6 

.0041  41.  E3   06  45.  E3  1.  47.  E3  *  1-3  PLASTICITY 


115*  K-l-1  FLOOR  SHELLS 

ie75  0.  2  *  K-2-1 
I  1  5  *  K-l-2  STIFFENER  SHELLS 
. 125  0   2  4  K-2-2 


1 

1 

2 

1 

2  0  4  tt-1-1 

0.  6.  0.  3. 

1  0 

410 


*  M-2A-1 
1  4  M-5-1 


4 

0 
0 

o. 

1  0 

410 
4  2 


4  M-2A-2 
1  *  M-5-; 


BOUNDS 


0 

0 

2 

0 

1 

2 

410 

4  b 


0   0.  0 
S  N-l-1 

2  4  $  P-l-1  BOUNDS  1 

Q-l-1 

R-l-1 

4  M-l-2 
i.      3.  9 

0   0.  0 

4  N-l-2 
2  2  6  4  P-l- 
4  Q-l-2 
$  R-l-2 

2  5  M-l-3 
6   0   . 9375  4  M-2A-3 

3  0.  3.   9375  4  M-3-3 
0  90   0.  0  1  4  M-5-3 

4  N-l-3 

2  6  4  P-l-3  BOUNDS  3 

Q-l-3 

R-l-3 

4  M-l-4 

3   3.  3594  4  M-2A-4 

0.  3  35^4   9375  4  M-3-' 

0   0   0  1SM-5-4 

s  N-l-4 

2  6  4  P-l-4  30UNDS  4 

Q-t-4 

R  -  1  -  4 

4  M-l-5 
6   2.  6406  3   4  M-2A-5 
110   2  6406   9375  4  M-3-: 
2  0  0   0   0  14  M-5-5 
410  4  N-l-5 

4  6  2  3    4  P-l-5  30UNDS  5 
0  i    Q-l-5 
0  J  R-l-5 


0  * 
0  4 
2  2 

0  i 

1  3 

2  0 
410 
4  3 
0  4 
0  $ 
2  2 
0 
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FLUID  MASS  DATA  FOR  TRANSVERSE 


TIFFENE? 


0  242  150 

0  0 

1.  152E-3 

59  E3^ 

T  T  F  T 

T  F  F  T 

F  F  T  F 

F  F  T  F 

F  T  F  F 

F  T 

TEE.  FLU  T 

EE.  GEO  TEE.  COR 

10 

1.   1.   1. 

1. 

MESH 

14  0  0  0 

17  8  2 

2  4  0  0  0 

2  9  9  3 

3  4  0  0  0 

3  9  10  4 

4  4  0  0  0 

4  10  115 

5  4  0  0  0 

5  11  12  6 

6  4  0  0  0 

7  13  14  a 

7  4  0  0  0 

8  14  15  9 

8  4  0  0  0 

9  15  16  10 

9  4  0  0  0 

10  16  17 

1  1 

10  4  0  0 

0 

11  17  18 

12 

114  0  0 

0 

13  19  20 

14 

12  4  0  0 

0 

14  20  21 

15 

13  4  0  0 

0 

15  21  22 

16 

14  4  0  0 

0 

16  22  22 

17 

15  4  0  0 

0 

17  23  24 

18 

16  4  0  0 

0 

19  2  5  26 

20 

17  4  0  0 

0 

20  26  27 

21 

ie  4  0  0 

0 

21  27  23 

19  4  0  0 

0 

22  23  29 

23 

20  4  0  0 

0 

23  29  30 

24 

21  4  0  0 

0 

25  31  22 

-1  L 
CU 

22  4  0  0 

0 

26  32  22 

27 

23  4  0  0 

0 

27  33  3  4 

23 

24  J  0  0 

0 
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29 

34  35 

29 

25 

4  0  0 

0 

29 

35  36 

30 

26 

4  0  0 

0 

31 

37  38 

22 

27 

4  0  0 

0 

32 

38  39 

33 

28 

4  0  0 

0 

33 

29  40 

34 

29 

4  0  0 

0 

34 

40  41 

35 

30 

4  0  0 

0 

35 

41  42 

36 

31 

4  0  0 

0 

37 

43  44 

38 

32 

4  0  0 

0 

38 

44  45 

39 

33 

4  0  0 

0 

39 

45  46 

40 

34 

4  0  0 

0 

40 

46  47 

41 

35 

4  0  0 

0 

41 

47  43 

42 

36 

4  0  0 

0 

43 

49  50 

44 

37 

4  0  0 

0 

44 

50  51 

45 

38 

4  0  0 

0 

45 

51  52 

46 

39 

4  0  0 

0 

46 

52  53 

47 

40 

4  0  0 

0 

47 

53  54 

48 

41 

4  0  0 

0 

49 

55  56 

50 

42 

4  0  0 

0 

50 

56  57 

51 

43 

4  0  0 

0 

51 

57  53 

en 

-  — 

44 

4  0  0 

0 

52 

58  59 

53 

45 

4  0  0 

0 

53 

59  60 

54 

46 

4  0  0 

0 

55 

61  62 

56 

47 

4  0  0 

0 

56 

62  63 

57 

48 

4  0  0 

0 

57 

o3  64 

58 

49 

4  0  0 

0 

53 

64  65 

59 

50 

4  0  0 

0 

59 

65  6o 

60 

51 

4  0  0 

0 

6  12  77  67 

52 

4  0  0 

0 

67 

77  73 

68 

53 

4  0  0 

0 

6a 

79  7n 

69 

54 

4  0  0 

0 

69 

~9    80 

70 
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55  4    0 

70  50 

56  40 

71  51 

57  4    0 

72  82 

58  4    0 

73  53 

59  4    0 

74  84 

60  4    0 

75  35 

61  4    0 
12    18 

62  4    0 

77  87 

63  4    0 

78  58 

64  4    0 

79  59 

65  4    0 

80  90 

66  4    0 

81  91 

67  4    0 

82  92 

68  4    0 

83  93 

69  4    0 

84  94 

70  4    0 


0  0 

81  71 

0  0 

tjcl  /  C 

0  0 

83  73 
0 


u 

84 
0 


74 
0 


35    75 
0    0 

86  76 
0    0 

87  77 
0    0 

88  78 
0    0 

89  79 
0    0 

90  SO 
0    0 


91 
0 

92 
0 


81 
0 

82 
0 


93  83 

0    0 

94  34 
0    0 

95  85 
0    0 


85    95    96    86 
71     4    0 


18    24 

72  4    0 

87  97 

73  4    0 

88  98 

74  4    0 

89  99 

75  4    0 

90  100 

76  4  0 
101 
4  0 
102 
4  0 
103 
4    0 


0 
87 


U 

97 

0    0 
98    88 

0    0 


91 
77 
92 
78 
93 
79 

94  104 

80  4    0 

95  105 

81  4    0 
24    30 

82  4    0 

97  107 
S3    4    0 

98  108 

84  4    0 

99  109 

85  4    0 


99  89 
0    0 

100  90 
0  0 
101 
0  0 
102 
0  0 
103 
0  0 
104 
0  0 
105 
0  0 
106  96 
0    0 

107    97 
0    0 

108  99 
0    0 

109  99 

0  0 

1  10 
0    0 


91 


92 


93 


94 


95 


100 


161 


100  1 1 0  111  101 

86  4  o  o  o- 

101  111  112  102 

87  4  0  0  0 

102  112  113  103 

aa  4  o  o  o 

103  113  114  104 

89  4  0  0  0 

104  114  115  105 

90  4  0  0  0 

105  115  116  106 

91  4  0  0  0 
30.  36  1  17  107 

92  4  0  0  0 

107  117  118  108 

93  4  0  0  0 

108  118  119  109 

94  4  0  0  0 

109  119  120  110 

95  4  0  0  0 

110  120  121  111 

96  4  0  0  0 

111  121  122  112 

97  4  0  0  0 

112  122  123  113 

98  4  0  0  0 

113  123  124  114 

99  4  0  0  0 

114  124  125  115 

100  4  0  0  0 

115  125  126  116 

101  4  0  0  0 
36  42  127  117 

102  4  0  0  0 
127  128  118 
4  0 

129  129  119 
4  0 


1  17 
103 

ue 

104 
119 

105 
120 
106 
121 
107 


108 
123 
109 
124 
110 
125 
1  1  1 
42  4B 
1  12  4 
127 
1  12 
125 
1  14 
129 
1  15 
130 


129 
4  0 
130 
4  0 
131 
4  0 
132 
4  0 
133 
4  0 
134 
4  0 
135 
4  0 

i; 
o 

137 
4  0 
133 
4  0 
139 
4  0 
1  40 


128 
0  0 
129 
0  0 
130 
0  0 
131 
0  0 
132 
0  0 
133 

0  0 

1  34 
0  0 
135 
0  0 
136 
0  G 

17  127 
0  0 
138 
0  0 
139 
0  0 
140 
0  0 

141 


120 


121 


122 


123 


124 


125 


126 


129 


130 


131 


162 


1 1  = 

4  0 

0  0 

121 

141 

142 

1  22 

1  17 

4  0 

0  0 

132 

142 

143 

123 

113 

4  0 

0  0 

132 

1-13 

144 

124 

119 

4  0 

0  0 

134 

144 

145 

125 

120 

4  0 

0  0 

1 35 

145 

146 

126 

121 

4  0 

0  0 

48  54  147  137 

122 

4  0 

0  0 

137 

147 

148 

123 

123 

4  0 

0  0 

123 

148 

149 

129 

124 

4  0 

0  0 

139 

149 

150 

140 

125 

4  0 

0  0 

140 

150 

151 

141 

126 

4  0 

0  0 

141 

151 

152 

142 

127 

4  0 

0  0 

142 

152 

153 

143 

123 

1  0 

0  0 

143 

153 

154 

144 

129 

4  0 

0  0 

144 

154 

155 

145 

130 

4  0 

0  0 

145 

155 

156 

146 

131 

4  0 

0  0 

54  60  157  147 

132 

4  0 

0  0 

147 

157 

153 

143 

133 

4  0 

0  0 

14S 

153 

1  59 

1  49 

134 

4  0 

0  0 

149 

159 

160 

150 

135 

4  0 

G  0 

150 

160 

161 

151 

136 

4  0 

0  0 

151 

161 

162 

152 

137 

4  0 

0  0 

152 

162 

163 

153 

133 

4  0 

0  0 

153 

163 

164 

154 

139 

4  0 

0  0 

154 

164 

165 

155 

140 

4  0 

0  0 

155 

165 

16o 

156 

141 

4  0 

0  0 

60  £ 

)6  I* 

,7  157 

142 

4  0 

0  0 

157 

167 

168 

153 

143 
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